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flux with support at infinity to local Calabi-Yau geometries in type IIB. This allows us to 
apply the formalism of Ooguri, Ookouchi, and Park [arXiv:0704.3613j to demonstrate that, 
by tuning the flux at infinity, we can stabilize the dynamical complex structure moduli in a 
metastable, supersymmetry-breaking configuration. Moreover, we argue that this setup can 
arise naturally as a limit of a larger Calabi-Yau which separates into two weakly interacting 
regions; the flux in one region leaks into the other, where it appears to be supported at infinity 
and induces the desired superpotential. In our endeavor to confirm this picture in cases with 
many 3-cycles, we also compute the CIV-DV prepotential for arbitrary number of cuts up to 
fifth order in the glueball fields. 
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1 Introduction 

Over the last years much progress has been made in studying flux compactifications in string 
theory; see [1,2] for recent reviews. By now there is strong evidence that there is a huge 
number of supersymmetric vacua with negative cosmological constant in which all scalar 
moduli are stabilized, the so called landscape of string theory. Typical constructions start 
with a warped Calabi-Yau compactification of type IIB string theory to four dimensions. Some 
of the scalar moduli are stabilized by the addition of fluxes through the compact cycles of the 
internal manifold and others by various quantum effects. Since supersymmetry is broken in 
the real world, to make contact with phenomenology it is necessary to extend the previous 
constructions to non-supersymmetric (meta) stable vacua with small positive cosmological 
constant. For this we need to understand the mechanism of supersymmetry breaking in 
string theory. By now several methods of supersymmetry breaking for string vacua have been 
proposed, such as the introduction of anti-branes [3,4], or simply the existence of metastable 
points of the flux-induced potential [5,6]. The main drawback of these constructions is that, 
in most cases, they are not under complete quantitative control. 

While the question of supersymmetry breaking should be ultimately understood in an hon- 
est compactification, that is in a theory including gravity in four dimensions, it is technically 
easier to study simpler systems where the gravitational dynamics has been decoupled from the 
gauge theory degrees of freedom. This typically happens in the limit where a local singularity 
develops in the Calabi-Yau manifold. In such a situation all the interesting dynamics related 



to the degrees of freedom of the singularity takes place at energy scales much lower than the 
four dimensional Planck scale. Assuming that supersymmetry breaking is related to these 
light degrees of freedom, it is then possible to zoom in towards the singularity and forget 
about the rest of the Calabi-Yau. This leads us to the study of supersymmetry breaking 
and string phenomenology in the context of local Calabi-Yau geometries possibly with the 
addition of probe D3-branes [7-1 1] . 

Meanwhile a new important aspect of supersymmetry breaking in gauge theories was de- 
veloped after the discovery of Intriligator, Seiberg and Shih [12] that even simple supersym- 
metric gauge theories can exhibit dynamical supersymmetry breaking in metastable vacua. 
From a phenomenological point of view this possibility is quite attractive, and a lot of activity 
has been concentrated around extensions of the ISS model and various related string theory 
constructions [13-22] (see also [23]). A certain class of gauge theories where supersymmetry 
breaking in metastable vacua can be studied with good control is that of M = 2 gauge theo- 
ries perturbed by a small superpotential, initiated by [24]. In such theories the exact Kahler 
metric on the moduli space is known, which allows one to compute the scalar potential pro- 
duced by the perturbation of the theory by a small superpotential exactly to first order in the 
perturbation. It was shown that generically there are metastable supersymmetry breaking 
vacua generated by appropriate superpotentials. We will refer to this as the OOP mechanism 
for supersymmetry breaking in M = 2 theories. 

String theory in a local Calabi-Yau singularity realizes geometric aspects of supersym- 
metric gauge theories. In particular the question of supersymmetry breaking in these two 
systems should be related. The first goal of our paper is to make this connection more pre- 
cise by giving a geometric realization of the OOP supersymmetry breaking mechanism in IIB 
on a local Calabi-Yau singularity. To realize OOP one first has to engineer the (IR of the) 
M = 2 gauge theory and then to find a way of introducing the appropriate superpotential. 
The first step is achieved by the standard geometric engineering of M = 2 gauge theories by 
IIB on noncompact Calabi-Yau manifolds [25,26]. It is well known that the moduli space of 
the Calabi-Yau compactification encodes the geometry of the Coulomb branch of the gauge 
theory and that the Seiberg- Witten solution can be rederived by the complex geometry of the 
Calabi-Yau. 

The introduction of superpotential to this system is less straightforward and, to our knowl- 
edge, has not been studied in the literature before, in this context. Our main proposal is that 
the superpotential can be introduced by turning on 3-form flux in the local Calabi-Yau, which 
is not piercing its compact 3-cycles, but which is growing in the noncompact direction of the 
Calabi-Yau. In other words, it is flux which has support at infinity. While this flux is not 



directly piercing the compact cycles we show that, once appropriately regularized, it does 
introduce an effective superpotential for the complex structure moduli, which is generaliza- 
tion of the usual Gukov-Vafa-Witten superpotential [27-29] to 3-form flux with noncompact 
support. This is a way to introduce a general superpotential in a geometrically engineered 
M = 2 gauge theory. In particular, we explain that in certain cases it is possible to engineer 
the OOP-type superpotential, which guarantees the existence of metastable, supersymmetry 
breaking vacua for the complex structure moduli. 

The second goal of our paper is to find a "natural" way to generate the supersymmetry 
breaking flux configurations described above, starting from a more standard setup. In this 
process, we also clarify the meaning of flux which has noncompact support and the various 
subtleties related to io The natural interpretation of the flux described in the previous 
paragraph emerges once we embed the previous supersymmetry-breaking local singularity into 
a bigger IIB compactification with standard flux of compact support. As shown in figure HJ 
the physical idea is to start with a Calabi-Yau manifold with a set of three-cycles which are 
isolated from the other three-cycles by a large distance. We turn on 3-form flux on all cycles 
except for the isolated set. While the flux that we have turned on is not piercing the isolated 
cycles, it does leak into their regional and it produces a potential for their complex structure 
moduli. In the limit where the distance between the two sets of cycles of the Calabi-Yau 
becomes very large, which we will refer to as the factorization limit, the flux leaking towards 
the isolated set will start to look like the flux coming from infinity, mentioned in the previous 
paragraph. In this way we manage to embed the scenario of the previous paragraph into a 
well defined system. While this factorization idea should work even in the case where the 
total Calabi-Yau is a compact manifold divided into two parts, U in this paper we will only 
analyze it in the local case, as it is technically easier. 

As a check of this, we consider the example of a local Calabi-Yau based on a hyperelliptic 
Riemann surface. In this case the factorization can be studied more explicitly. Matrix model 
techniques can be used to compute the prepotential in the factorization limit. Our results 
verify the general intuition of the last paragraph. 

The plan of this paper is as follows. In section [2J we review some general aspects of 
IIB flux compactifications and the potentials that can be generated by noncompact fluxes in 
the local limit. We also discuss the general mechanism by which such fluxes may be used 



For example the fact that naively the total energy density in four dimensions diverges. 

2 This means that the 3-form field strength is nonzero in the region around the isolated set of 3-cycles, but 
once integrated over one of these 3-cycles the integral is zero. 

3 Because of no-go theorems [30-32], in such compact setups one will need extra ingredients such as O- 
plancs, which we do not consider in the present paper. 
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Figure 1: Idea of the paper: In I we start with a generic Calabi-Yau with flux 
piercing through some of its 3-cycles, while making the distance between the cycles 
with and without flux very large in II. This is seen as flux from infinity in the left 
sector without compact flux in III, and generates an OOP-like potential in that 
sector. 



to stabilize complex structure moduli in metastable supersymmetry-breaking configurations. 
After this, we turn in section [3] to the study of local Calabi-Yau geometries based on Riemann 
surfaces, providing a more detailed description of the generation of metastable vacua in this 
context and providing an explicit example. Section H] then addresses the second point of this 
paper, namely the ability to obtain our local geometries with noncompact fluxes from larger 
Calabi-Yau with compact ones by taking a factorization limit. Section [5] supplements this 
general discussion by providing an explicit demonstration, using matrix model techniques, of 
the factorization limit in the class of local geometries studied in section [31 Finally, we finish 
in section [6] with some concluding remarks concerning the generalization of our story to other 



M = 2 contexts, such as M and F-theory compactifications. Some supplementary material 
and technical details are contained in four appendices. In appendix O, in particular, we study 
prepotential for the Cachazo-Intriligator-Vafa/Dijkgraaf-Vafa geometry using matrix model 
as a step towards confirming the above scenario of realizing metastable vacua by factorization. 
We compute the prepotential for an arbitrary number of cuts up to fifth order in glueball field. 

As this paper was being prepared for publication, a paper appeared [33] where a similar 
system was studied as an example from a different perspective. 

2 IIB Compactifications with Flux at Infinity 

2.1 Compact Calabi-Yau 

Compactification of IIB on a Calabi-Yau threefold M. leads to an M = 2 supergravity theory in 
4d. The number of vector multiplets is h 2,1 , their scalar components correspond to the complex 
structure moduli of A4. We also have h 1 ' 1 + 1 hypermultiplets, whose scalars correspond to 
the Kahler moduli of M. and the axion-dilaton. The two sets of multiplets are decoupled, and 
in the rest of the paper we will concentrate on the dynamics of the vector multiplets. 

A Calabi-Yau threefold has a nowhere vanishing holomorphic (3, 0) form Q which is unique 
up to scale. Consider a symplectic basis of 3-cycles {A 1 , Bj} with I, J = 0, 1, ... , h 2 ' 1 . We 
define the periods of Q as 

x 1 = I n, Fi= I a (2.i) 

J A 1 JBi 

The ^4-periods X 1 are projective coordinates on the complex structure moduli space of Ai, 
and the Fi are functions of X 1 . The metric on the complex structure moduli space is special 
Kahler and the Kahler potential is given by 

K = -log(i [ttAtt\ . (2.2) 

This is an exact result which does not receive any a' or g s corrections. 

The easiest way to lift the phenomenologically unrealistic moduli space of such compacti- 
fications is to turn on fluxes through the compact cycles of the Calabi-Yau. In the case of IIB 
we can turn on RR and NS-NS 3-form flux F% and H3 through the 3-cycles of the threefold. 
This generates a superpotential for the complex structure moduli [27-29] given by 

W= /GsAfi, (2.3) 



where G 3 = F 3 — tH 3 and t = C$ + i/g s - The scalar potential is computed by the standard 
Af = 1 supergravity expressions 

V = e k (G al D a WDbW - 3\W\ 2 ) , (2.4) 

where G a i is the metric on the moduli space derived from the Kahler potential K, and where 
we have introduced the Kahler covariant derivative D a W = d a W + (d a K)W. 

The F 3 and H 3 fluxes generate charge for the F$ form via a Chern-Simons coupling in the 
lOd IIB supergravity action. The F 5 flux has nowhere to end, so we are lead to the tadpole 
cancellation condition for IIB compactifications 

^Jf 3 AH 3 + Q D3 = } (2.5) 

where Qr> 3 receives positive contribution from probe D3 branes and negative contribution 
from induced charge on D7 and orientifold planes. 

2.2 Local Limit 

It is well known that there are points on the complex structure moduli space of Calabi- 
Yau manifolds where the manifold develops a singularity [34]. The simplest example is the 
conifold singularity, where we have a 3-cycle whose size goes to zero. More generally, a more 
complicated set of cycles may become very small in some region of the moduli space. As 
we approach this region, the local dynamics of the singularity decouples from the rest of the 
fields. What this means is that in 4d the typical energy scale for the dynamics of the fields 
corresponding to the singularity becomes much smaller than any other scale, in particular 
much smaller than the Planck mass M p in 4 dimensions. In this sense, the dynamics of 
the singularity is decoupled from gravity. Moreover to study the relevant dynamics, we can 
zoom in close to the singularity and forget about the rest of the Calabi-Yau. In this limit 
the Calabi-Yau looks noncompact, and it becomes technically easier to study the low energy 
dynamics. 

A typical example of such a local Calabi-Yau is a complex manifold of the form of a 
hyper surf ace in C 4 

M : uv-F{x,y) = 0, (2.6) 

where F(x, y) is a polynomial. In this case the holomorphic 3-form is 

du 
Q = — AdxAdy. (2.7) 

u 



In this expression the indices a, b run over complex structure moduli, Kahler moduli and the axion-dilaton. 
We denote by K the total Kahler potential for all moduli and by K, as in (|2.2[) . the one for the complex 
structure moduli alone. 



By taking the local limit to go from a compact Calabi-Yau to a noncompact one, the structure 
of special geometry described above reduces to what is called rigid special geometry [35,36], 
which is relevant for the low energy dynamics of M = 2 gauge theories. In this case the 
Kahler potential reduces to 

K = i J QAQ, (2.8) 



and the Kahler covariant derivative Di reduces to the ordinary derivative <9j. 

An important point is the distinction between normalizable and non-normalizable complex 
structure moduli in the case of noncompact Calabi-Yau manifolds. To be more precise let us 
consider the example ( 12. 6j) . The coefficients {f} of the polynomial P(x,y) characterize the 
complex structure of A4, so {£'} are the complex structure moduli of M.. However not all 
of them are dynamical. Some of them control the complex structure of 3-cycles which are 
localized in the "interior" of the singularity and are dynamical, while others describe how 
the singularity is embedded in the bigger Calabi-Yau and become frozen when we take the 
decoupling limit. 

To determine if a specific complex structure modulus t is dynamical or not, one has to 
compute the corresponding Kahler metric 

g tI = d t dlK = id t W t j O A H. (2.9) 

If this expression is finite, then the modulus t is dynamical, otherwise it is decoupled and 
should be treated as a parameter of the theory. We will refer to the first set of moduli as 
normalizable and to the second as non-normalizable. 

2.3 Adding Flux 

As in the compact case, the addition of fluxes to the local Calabi-Yau introduces a superpo- 
tential for the moduli. The dynamics of the Kahler moduli and the dilaton decouple, and we 
can concentrate on the normalizable complex structure moduli. The superpotential is still 
given by (12.31) . but now the scalar potential is computed by the rigid M = 2 expression 

V = G i3 d i WtyW= f G 3 A*G~ 3 . (2.10) 

Since we are in a noncompact Calabi-Yau it is not necessary to impose the tadpole cancellation 
condition. Instead, the quantity 

[f 3 AH 3 (2.11) 

represents the F5 flux going off to infinity and remains constant as we vary the moduli. We 
will use this to simplify the potential in the next section. 
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In most treatments of fluxes in noncompact Calabi-Yau manifolds the assumption is made 
that the flux is threading the compact cycles of the singularity and is going to zero at infinity. 
As we explained in the introduction the goal of our paper is to study the dynamics in the case 
where the flux is actually coming in from infinity and is not supported on the compact three- 
cycles. Of course, in a local singularity inside a bigger compact Calabi-Yau, what is meant 
by infinity is the rest of the Calabi-Yau and we should think of flux coming from infinity as 
flux leaking towards the singularity from the other compact cycles. 

More precisely, in a noncompact Calabi-Yau we consider the vector space H 3 (Ai) of 
harmonic 3- forms which do not necessarily have compact support, so they can grow at infinity. 
The harmonic 3- forms of compact support form a linear subspace H^ pct (Ai ) C H 3 (At ) . There 
is a natural way to define the complement subspace H^ (Ai ) C H 3 (Ai ) as the harmonic forms 
with vanishing integrals on the compact 3-cycled_|. Then we have the decomposition 

H 3 (M) = H^M) © H 3 cva {M). (2.12) 

We will also refer to the forms in H 3 t (Ai) as harmonic 3-forms with compact support and 
to those in H^Ai) as 3-forms with support at infinity. 

Now we want to consider the case where the 3-form field strength that we have turned on 
has support at infinity 

G 3 e Hl(M), (2.13) 

which means that G 3 has zero flux through the compact cycles 



G 3 = / G 3 = 0. (2.14) 

A}- JBi 

The intuitive picture that one should keep in mind, is that this flux at infinity represents 
usual flux piercing other 3-cycles which are very far away from the singularity in the big 
Calabi-Yau. As we will see in more detail in the next section, in this case and if one zooms 
into the local singularity it is a good approximation to treat the flux from the distant 3-cycles 
as flux which "diverges" at infinity. In other words both H^Ai) and H 3 pct (Ai) correspond 
to the usual H 3 pct (Ai) of the bigger Calabi-Yau M. in which the singularity A\ develops. 

What is maybe more surprising is that the 3-form flux G 3 with support at infinity generates 
a potential for the complex structure moduli of the singularity At, even though it is not 



5 We should clarify that we arc not interested in the most general harmonic 3-form with noncompact 
support, but only in a restricted subset characterized by 3-forms which grow in a "controlled" way at infinity. 
This means that we want to consider forms which have at most a "pole" of finite order at infinity, and not 
essential singularities. This statement has a nice interpretation in the example where we have a local Calabi- 
Yau based on a Riemann surface that we will study later. Another way to state this restriction is that we 
will consider harmonic 3-forms on a local Calabi-Yau which do have a lift to the original Calabi-Yau that we 
started with before we took the local limit near its singularity. 

9 



directly piercing the compact cycles of Ai, as can be seen from (12.141) . Our starting point for 
the computation of this potential is the energy stored in the 3-form field 



V = G 3 A*G 3 . (2.15) 

Since G 3 has noncompact support, this is a divergent integral meaning that the energy of the 
flux is infinite. This was to be expected and is not really a problem, since we are interested 
in the changes of this energy as we vary the sizes of the 3-cycles in the neighborhood of the 
singularity. We would like to throw away the divergent, moduli independent piece of this 
quantity and keep the finite, moduli dependent one. A nice way to achieve this is to use the 
fact that the net F 5 form flux leaking off at infinity, being a topological quantity, has to be 
kept constant as we vary the moduli. It is easy to show that we can write 

G 3 AG^=(r-r) [f 3 AH 3 , (2.16) 



and the left hand side must be constant for the reason we explained. Since it is a constant 
we can subtract it from the potential and define 



V = J G 3 A *G 3 - / G 3 A G 3 . (2.17) 

It is easy to show that this is equal to 

V = jG 3 A*'Gj, (2.18) 

where G 3 is the imaginary anti-self dual part of the G3 flux 

*G 3 =-iG 3 . (2.19) 

The expression (" 12. 18[) is the finite and moduli dependent piece of the potential (12.151) . 

2.4 Simplifying the Potential 

In this section we simplify the expression (12.181) for the potential. In general we have the 
following relation between the Hodge decomposition and the * operator on a threefold 

* H 3 >° = -zH 3 * , *H^ 2 = -iH^\ 

(2.20) 

Before we proceed we would like to analyze the relation between the decomposition (12.121) 
and the Hodge decomposition. In general we have the following decomposition.] 

H 3 (M) = H 3 J> © H 3 cp ° ct © H 2 J © H 2 p \ t © {c.c.}. (2.21) 



6 Again, we are only considering a certain subset of all harmonic 3-forms with noncompact support, as 
explained in footnote [5] 
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Harmonic forms in H^ ct have compact support, w 



rile those in H™ do not, and are chosen 



a 



to have vanishing ^.-periods on the compact cycled. Since we do not want to break super- 
symmetry explicitly by the boundary conditions of the system, we want our configuration to 
be supersymmetric at infinity, which means that the flux at infinity has to be imaginary self 
dual so 

G 3 EH 2 J®H 2 J ct ®H 1 cp 2 ct . (2.22) 

where the subscript oo means that we have to consider the elements of the cohomology with 
noncompact support. We pick a basis 



E m G H^ , tti G H' t (2.23) 



with the following periods 



H m = 0, / ^ = 5}, 

A' J A* 

5 = K / O • = r • 

B, JBi 



(2.24) 



where r^- is the period matrix of the Calabi-Yau, and K im are holomorphic functions of the 
normalizable-complex structure moduli. 
The flux has an expansion of the form 

G 3 = T m E m + 11% + FH" (2.25) 

The parameters T m are fixed by the boundary conditions and have to be kept constant as we 
vary the normalizable moduli. We have also assumed that 



which means 



G 3 = G 3 = 0. (2.26) 

A i JBi 



T m E m + h? n d + V / % = 

J f J f _ J f_ (2.27) 

T m / E m + h? \ toj + V / n d = 0. 

JBi JBi JB % 



The first equation of (I2.27P implies that 



V = -h j . (2.28) 



7 A harmonic (p, <7)-form cannot have vanishing periods on all compact cycles unless it is identically zero. 
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and the second 

h '--k{^) 1(K "" Tm) - (229) 

As we explained before, only the imaginary anti-self dual part of the flux Gj = l l Vti con- 
tributes to the regularized potential and we have 

V = [ G 3 A G3" 



= -(K im T™)^—) (K jn T n ). (2.30) 

In this final expression the period matrix r y and Ki m are functions of the normalizable 
complex structure moduli, while T m 's have to be considered as constants which play the role 
of external parameters. This potential is in general very complicated and can have local 
nonsupersymmetric minima for appropriate choices of the parameters T m as we will explain 
latec. 

2.5 Properties of the Potential 

The potential (12.301) should look somewhat familiar as it shares the same basic structure as the 
scalar potential that arises when one adds a small superpotential to Seiberg-Witten theory. 
This connection can be made even more transparent by noting that Ki m can in general be 
written as a total derivative with respect to the special coordinates X 1 )q 

K im = ^.K m (X j ), X 1 = I Q. (2.31) 

With this notation, (12.301) takes the standard form 



where 

W eS (X k ) = T m K m (X k ) (2.33) 

is in fact proportional to the Gukov-Vafa-Witten superpotential induced by the flux G3. 



8 Although we do not discuss this in the present paper, from the viewpoint of flux compactification it is 
a natural generalization to consider fluxes through the compact 3-cycles, relaxing the condition (|2.26[) . Such 
flux will make additional contribution to the superpotential of the form N r Fi — aiX 1 , cti = J„ fl, which 
cannot be controlled by external parameters and makes realization of OOP-like vacua more difficult. 

9 One quick way to see this is to use the identity J M S m A difi — to derive Ki m ~ J QM A m A diQ, for 
a 2-form A m satisfying dA m — S m on the boundary (at infinity) of M. . Because the divergent contributions 
to A m at infinity can be chosen independent of the dynamical moduli, we can pull the derivative outside of 
everything. 

12 



2.5.1 OOP Mechanism 

Equation (j2.32p makes manifest the relation between our flux-induced potential (12.301) and 
that which arises in deformed Seiberg-Witten theory and allows us to utilize the technology de- 
veloped by Ooguri, Ookouchi, and Park [24] in that context] 10 ! for engineering supersymmetry- 
breaking vacua. In particular, if we want to realize a nonsupersymmetric minimum at some 
point P' ' in the moduli space, the OOP procedure tells us to first construct Kahler normal 
coordinates [38-40], around X 1 ^ 

oo 1 

z< = AX 1 + g i3 J2 ~ d h ■ ■ ■ d in T jhl2 AX h AX^ . . . AX 4 ", (2.34) 

n=2 

where AX 1 = X % — X 1 ^ and ~ means evaluation at X = X^°\ We then build the potential 
V in ( 12.321) from a superpotential W e s consisting of a linear combination of the z % 

W e ff = kiz\ kf. constant. (2.35) 

Stability can then be demonstrated by expanding V near p 

V = kihjgv + kikjR iJ kI z k z ! +0{z 3 ). (2.36) 

The curvature of special Kahler manifolds, of which the complex structure moduli space is an 
example, is positive definite at generic points. As a result, any potential of the form (12.321) 



that agrees with ( 12.351) near X % ^' to cubic order will engineer a nontrivial vacuum at X % ^ o 
One can obtain a nice physical picture for this mechanism by noting, as in [41], that the 
series (12.341) can be summed exactly and inserted into (I2.35P to yield 

W eS ~ e.X' + rriFi, (2.37) 

where e, and m? satisfy 

ei + m j f ji = 0. (2.38) 

From this, we see that the superpotential (12.371) built from Kahler normal coordinates is of 
precisely the form that we would have obtained had we instead simply turned on compactly- 
supported fluxes m % and e^ threading the cycles A 1 and £>;, respectively. The condition (I2.38p . 
however, combined with the requirement that Im r be positive definite, implies that e; and m l 



10 See also related work by Pastras [37]. 

11 For non-generic X^°\ the curvature may have a zero eigenvalue in which case higher order agreement 
with (12.35|) is required (that X % (°' is a stable for superpotential exactly equivalent to (|2.35j) will follow from 
the discussion below). 
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can never satisfy the condition e, + mPfji = that is required for preservation of the manifest 
M = 1 supersymmetry. 

It is well-known that the flux-induced superpotential (12.371) only breaks the full Af = 2 
supersymmetry spontaneously, though, so there is a second J\f = 1 in the game that is not 
manifest in this formalism. The relation f!2.38|) is, in fact, nothing other than the condition 
that the vacuum at X 1 ^ preserves precisely these non-manifest supersymmetries [41,42]. As 
such, the vacuum at X 1 (°> in the presence of the superpotential (12.351) is stable for a good 
reason — it is secretly supersymmetric! 

In general, our noncompactly supported fluxes will not generate potentials with W e s ex- 
actly equivalent to (12.351) . Rather, the W e fj's that arise are globally well-defined functions 
on the moduli spaccj which we can then tune to agree with ( 12. 35ft to cubic order within a 
neighborhood of the point X l( -°\ The delicate manner by which the superpotential (12.351) 
managed to realize a non-manifest M = 2 supersymmetry is crucially dependent on the full 
infinite series expansion about X 1 ^ so, by failing to exactly reproduce (I2.35p . we are able 
to explicitly break, at the level of the Lagrangian, the half of supersymmetry which would 
otherwise have been preserved by the vacuum at X t( ^\ Stability of the X 1 ^ 1 vacuum, on 
the other hand, depends only on the local behavior of W e g so our procedure will retain this 
property, leaving us with a locally stable supersymmetry-breaking vacuum. 

In the end, what we are doing to engineer a supersymmetry-breaking vacuum at X 1 ^ is 
actually a quite intuitive procedure. We first turn on a collection of noncompactly supported 
fluxes which explicitly break half of the M = 2 supersymmetry. We then tune these fluxes so 
that, near X % ^\ their interactions with the dynamical complex structure moduli mimic those 
of the compactly supported fluxes that would generate a vacuum at X 1 (°) which preserves the 
opposite half of supersymmetries. 

2.5.2 Supersymmetric Vacua 

In addition to possessing supersymmetry-breaking vacua when the T m are suitably tuned, the 
potential ( 12.32ft also typically contains a wide variety of supersymmetric vacua. As discussed 
in [24], these vacua fall into two different classes. First, because there is no flux directly 
threading the compact cycles, the energy cost associated with shrinking them is necessarily 
finite. Because the period matrix r^ diverges, the potential vanishes at these singular points 
and we obtain stable vacua which are in fact supersymmetric. 

This argument is of course rather crude because we are neglecting the new light degrees of 
freedom that enter as 3-cycles degenerate but, as is well-known, this is easily fixed. In particu- 



12 Contrast this with (|2.3T[) . which manifestly suffer from monodromies for constant (non-transforming) hi. 
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lar, the light D3 branes which wrap the degenerating cycles give rise to hypermultiplets [43,44] 
comprised of pairs of N = 1 chiral superfields Qi and Qi with bilinear superpotential couplings 
to special coordinates. For the simple case of degenerating A 1 cycles, the superpotential takes 
the schematic form 

W = W eS (X l ) + (QQ)iX\ (2.39) 

and allows a supersymmetric vacuum at X 1 = through condensation of QQ 

(QQ)i = -^? ( XJ = °) • ( 2 - 4 °) 

The second class of supersymmetric vacua correspond to solutions of the F-term equations 

diWrtiXi) = 0. (2.41) 

In general, there may be many solutions to these equations, as we will explore later in the 
example of section 13.21 

2.5.3 Lifetime of Supersymmetry-Breaking Vacua 

Because we have managed to achieve supersymmetry-breaking vacua while freezing all non- 
normalizable moduli, the energies Vo will in general be finite and independent of the cutoff 
scale Ao that we use to regulate the local geometry. This means that our vacua are truly 
metastable, even within this local model, and can decay to any of the supersymmetric vacua 
that exist in these models. Because the number the supersymmetric vacua is potentially large 
and their properties quite model-dependent, it is difficult to make general statements about 
the lifetime of our OOP vacua. Nevertheless, we recall here one observation from [24], namely 
that the decay rates will in general scale like 

e~ S , S ~ ^, (2.42) 

where Az is the distance in field space between the initial and final vacuum state and V + 
is the difference in their energies. By simultaneously scaling all T m by a common factor, 
rpm _^ e ym^ we can re tain our supersymmetry-breaking vacua while decreasing V + by the 
same factor, V+ — > eV + . In this manner, we see that, just as with OOP vacua in deformed 
Seiberg-Witten theory, these OOP flux vacua can be made arbitrarily long-lived 13 !. 



13 Because we should really think of the local Calabi-Yau as sitting inside some larger compact geometry, 
one important caveat to this statement of longevity is that the noncompact fluxes T m in reality derive from 
a suitable set of compact fluxes in the full Calabi-Yau. This means that there will be a series of quantization 
conditions that must be imposed that may affect the degree to which they may be tuned. 
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3 Metastable Flux Vacua in Local Calabi-Yau 

In the previous section, we saw that, starting from a compact Calabi-Yau and taking a 
decoupling limit, one ends up with a local Calabi-Yau with noncompact flux with support at 
infinity, which is nothing but the flux leaking from the rest of the full Calabi-Yau that have 
been decoupled, towards "our" local Calabi-Yau. Furthermore, this noncompact flux induces 
potential (12.301) for the complex structure moduli in the local Calabi-Yau. Depending on 
the noncompact flux, this potential can be very complicated and create nonsupersymmetric 
metastable vacua in the local Calabi-Yau; the OOP mechanism [24] reviewed in 12.5.11 tells 
us exactly how this can be done. In this section, we will take specific examples of local 
Calabi-Yau and demonstrate that one can generate such OOP vacua as IIB flux geometries. 



In subsection 13. 1[ we review constructions of typical local Calabi-Yau geometries, taking 
Seiberg-Witten and Dijkgraaf-Vafa geometries as examples. The focus will be on the form 
of the potential for the moduli which is induced by flux at infinity. We also make remarks 



on the gauge theory interpretation of the physics of these geometries. In subsection 13. 2\ we 
proceed to an explicit construction of metastable flux vacua in a Dijkgraaf-Vafa geometry, by 



tuning superpotential appropriately. In subsection 13.31 we estimate how much control of flux 
at infinity is required to create OOP vacua. 

3.1 Local Calabi-Yau Based on Riemann Surface 

A large group of examples of noncompact Calabi-Yau manifold in IIB is defined by an equation 
of the form 

uv-F(x,y) = 0, (3.1) 

where x,y can both be variables in C or C*. Compactifying on such a Calabi-Yau leaves 
M = 2 supersymmetry unbroken in four dimensions. Important roles in these Calabi-Yau's 
are played by the underlying one-dimensional complex curve in the x, y-plane defined by 
F(x,y) = [25,45]. In most of our examples this curve is smooth, and we will refer to it as 
the Riemann surface S. The total Calabi-Yau space will be named .Ms- The holomorphic 
3-form of A^s is given, e.g. for x, y € C, by 

^ du A dx A dv du , 

Q = xtix =—AdxAdy. 3.2 

Notice that the total threefold can be described as a local (or decompactified) elliptic 
fibration over the x, ?/-plane. Over generic points in the base x, y-space, its fiber is described 
by a hyperboloid satisfying the equation uv = fi where /i is nonzero, which may be viewed 
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as a decompactified compact torus, making its B-cycle very large. On the other hand, when 
(x,y) G E, the noncompact fiber degenerates into a cone uv = 0, which one obtains by 
decompactifying a pinched torus, corresponding to an A geometry. 

Many important properties of the noncompact Calabi-Yau threefold Ai^ have an inter- 
pretation in terms of the underlying Riemann surface E. For example, the compact 3-cycles 
{A 1 , Bj} in A^e are lifts of compact 1-cycles on E, which we denote by {A\ Bj}. If (x, y) G C 2 , 
these 3-cycles may be constructed by filling in a disk D in C 2 whose boundary dD is the 1- 
cycle on E. Now consider an S^-fibration over D where S 1 is the compact circle in the uv -fiber. 
Since this circle shrinks over E, the total 3-cycle has the topology of an S 3 . If one of the 
variables x or y is C*-valued, the disk D will be punctured. In such a situation differences 
of 1-cycles have to be considered. We will see an example of this shortly. Notice that the 
one-to-one correspondence between 3- and 1-cycles shows an equivalence between the complex 
structure moduli on .Ms and E. 

A basis of (2,l)-forms with compact support on A^s is given by derivatives of Q with 
respect to the normalizable complex structure moduli: {Qi = diQ}. If M.^, were compact, 
these derivatives <9j would be Kahler covariant derivatives Di on the moduli space. Being 
noncompact instead, the moduli space is described by rigid special geometry and, as we saw 
before, the covariant derivatives simplify into partial derivatives. Another reduction over the 
compact 3-cycles in the Calabi-Yau shows that all these compactly supported (2, l)-forms fij 
reduce to a basis of holomorphic 1-forms Ui on E. Similarly, (l,2)-forms <9jf2 in Ai-% reduce 
to antiholomorphic forms uJj on E. Ui satisfy the following relations, which are reductions of 
(E21: 

h\^ = ^ h\ B ^ = T -- (3 - 3) 

The relation between the 3-cycles/3- forms on A4s and the 1-cycles/ 1-forms on E through 
the trivial wf-fibration being understood, we can rewrite the various relations in section [2] in 
terms of the Riemann surface E. First of all, the holomorphic 3-form Q of A^s, which is given 
e.g. for x, y G C by ( 13. 21) . is easily seen to reduce to a meromorphic 1-form A = ydx on the 
Riemann surface in this case [25,45]. The special coordinates ( 12. ip parametrizing complex 
structure moduli are 

X 1 = ^-. f A, F t = -L / A, (3.4) 



2ni J Al ' 2iri j B . 



and the Kahler potential (12. 8p is given by 

K = i / AAA. (3.5) 



x; 
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Recall that, in the special coordinates {X 4 }, the moduli space metric takes a particularly 
simple form: 

/ d 2 K \ 

ds 2 = = dX l dXi = (Imr) w dX l dXi (3.6) 

KdX'dXo) K >3 K ' 

as can be shown using <9jA = Ui and the Riemann bilinear relation. 

Now we want to consider a very small deformation of the system breaking supersymmetry 
to M = 1, thus generating a potential V for the moduli. As we saw before, this can be 
accomplished by turning on 3-form flux G3 with support at infinity in the local Calabi-Yau. 
This flux an be thought of as leaking from the other part of the full compact Calabi-Yau, 
which has been frozen in the decoupling limit. We assume that the decoupling limit was taken 
consistently with the elliptic fibration structure; namely, we assume that the noncompact flux 
is supported at the asymptotic infinities of E, while being compact in the direction of the 
■uv-fibers. 

The basis of (2,l)-forms with noncompact support, {H m }, in the Calabi-Yau Ms descend 
to meromorphic 1-forms {£ m } on the Riemann surface E, satisfying the relations 

/ U = 0, f U = K im , (3.7) 

which are reductions of (I2.24p . Therefore, the 3-form flux with noncompact support, G3, on 
.Ms as given in ( 12.251) descends to a harmonic 1-form flux 

9 = 9H + gA, , N 

_ __ (3.8) 

g H = T m i m + h u t , g A = l' l uJi, 

which will have poles at the punctures (or asymptotic legs) of E. A 3-form flux G3 in .Ms 
induces superpotential (12.31) . which reduces to an integral on E: 

W= [ gAX, (3.9) 

while the associated scalar potential (12. 18[) reduces to an integral on E: 



V = g A AgI. (3.10) 

If we require the condition (12.261) that the flux (13. 8p is zero through compact 3-cycles of 
A^sW) which translates into 

r g= f g = 0, (3.11) 

A* JBi 



then by the exactly same argument we did for general Calabi-Yau's in the previous section 
now reduced to the Riemann surface £ (or simply by borrowing the result (J2.30P ). we can 
rewrite (13. IQj) in terms of periods on S: 



V = - A {K im T^)\±-XK^. 



l j™ j 



(3.12) 



For convenience, the relation between 3- and 1-forms in A^s and £ is summarized in Table 

m 





special forms 


noncompact rlux 
inducing superpotential W 


compact flux 
entering potential V 


Mx 


Q e H 3 '°(Ms) 


T m E m e H 2 J(Mx) 


Gt=m i eHl£ t {Mx) 


£ 


A G H^CE) 


T m U g if^°(E) 


gJ = J l uT i e fl^t(S) 



Table 1: The summary of the relation between 3-formd_| in Calabi-Yau A^s and 
1-forms on Riemann surface E 



Now we will turn to more specific examples of local Calabi-Yau geometries based on 
Riemann surfaces which have been studied in the context of string theory. 

3.1.1 Seiberg-Witten Geometries 

An illustrative example of the general Calabi-Yau's (13.11) is given by SU(N) Seiberg-Witten 
(SW) geometries. In type IIB, these correspond to compactifications on noncompact Calabi- 
Yau threefold A^sw defined by 



Afsw : uv - F sw (x, y) = 0, x G C, y G C, 



(3.13) 



where the underlying Riemann surface Esw is a hyperelliptic curve 



^sw 



1 



F sw (x,2/) = A 7V [y + ~] -Pn(x) = 



and Pn(x) is a polynomial of degree N with the coefficient of x being zero: 



N 



N 



Pn(x) = Y\(x -act), ^ a i = 0- 



i=i 



i=l 



(3.14) 



(3.15) 



15 As explained in footnote El we do not mean here that S m 's span a complete basis of 3-forms in Air, with 
support at infinity; we are only considering a certain subset of all 3-forms diverging at infinity, which are the 
lifts of meromorphic 1-forms £ m on M. . 



19 



The coefficients of P^(x), or equivalently aij, are normalizable moduli, while A is a fixed 
parameter. The holomorphic 3- form on A^sw is ^sw = — A dx A — and reduces to [25] 

Asw = x — (3.16) 

y 

on the Riemann surface £sw- 

Type IIB string theory compactified on the Calabi-Yau (13. 13[) without flux geometri- 
cally engineers [25,46] an Af = 2 Seiberg-Witten theory [35,47]. In particular, the SU(N) 
Seiberg-Witten curve of gauge theory [48,49] is geometrically identified with the curve (13.141) 
underlying the Calabi-Yau. A T-duality along the compact circle in the uv-ftoer, followed by 
a lift to M-theory, translates [50] this geometry into a system of an M5-brane which wraps 
the Riemann surface Esw and fills IR 3,1 . In the IIA limit, this system is related to a Hanany- 
Witten type brane configuration in type IIA, where one has two NS5-branes with A" D4-branes 
stretching between them [51,52]. From this last IIA/M-theory point of view, it is easy to see 
the relation of the system to M = 2 SU(N) super Yang-Mills as the worldvolume theory on 
the D4-branes. In particular, a;'s correspond to the eigenvalues of the adjoint scalar $ on 
the Coulomb branch. In passing we also note that the geometries (13. 13[) are related to toric 
geometries in IIA by mirror symmetry [26,53]. 

Now let us look at the homological structure of the Seiberg-Witten geometry ( 13.131) . fo- 
cusing on the relation between 1-cycles on the hyperelliptic curve Ssw H3.14J) and 3-cycles in 
the Calabi-Yau A^sw (13. 131) . The Riemann surface S sw may be compactified by adding two 
points at infinity. If we represent the curve (13.141) as a two-sheeted x-plane branched over 
2N points, those infinities correspond to x = oo on the two sheets. It is thus a hyperelliptic 
curve of genus A^ — 1 with two punctures. Therefore, its first homology ifi(Ssw) is formed by 
N — 1 pairs of compact A and 5-cycles, (A 1 , Bj), i,j = 1, . . . , N — 1, with in addition a closed 
1-cycle A°° around one of the punctures which is dual to an open 1-cycle B^ connecting the 
two points. How can these 1-cycles be lifted to 3-cycles in A^sw? The fact that t/eC* means 
that A-cycles on Egw are not contractible on the x, y-plane (recall that we are regarding E$w 
as embedded in the x , ?/-plane) . Instead, compact A-cycles in the noncompact Calabi-Yau 
threefold will reduce to differences of A-cycles on £sw- Indeed, notice that a point on the 
1-cycle A 1 and one on another 1-cycle — A\ with opposite orientation, are connected by a P 1 
in the Calabi-Yau. The resulting 3-cycle therefore has the topology of S 2 x S 1 . For the B- 
cycles this subtlety does not arise, and compact B-cycles in the Calabi-Yau have S 3 topology 
and reduce to compact 1-cycles connecting the two hyperelliptic planes. See Figure [2j This 
discussion is equivalent to page 10 of [25], and in particular shows the equivalence between 
3-cycles on the Calabi-Yau's and 1-cycles on the Seiberg-Witten curve. 
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.A-cycle 




B-cycle 




Figure 2: The relation between 3-cycles in the Calabi-Yau .Msw and 1-cycles on 
the Riemann surface Esw for the Seiberg-Witten geometry. For the .4-cycle, by 
fibering S 1 over the line segment whose endpoints are at a point on A 1 and a point 
on — A\ one obtains S 2 . By moving the endpoints over A 1 and — A\ one obtains 
S 2 x S . For the £>-cycle, similarly moving the S 2 ending on Bi, one obtains S 3 . 



As seen in section 13.11 the complex structure moduli space is conveniently parametrized 
by the special coordinates ( 13.41) . which in the Seiberg-Witten case is conventionally denoted 
by a\ i = l,...,JV-l@ 

dy 



A 



1 



sw 



X- 



2m J a% 2m' J A i y 

As in (13.61) . the moduli space metric takes the special form for these: 

ds 2 = {\n\Tij)da % daK 



(3.17) 



(3.18) 



Using a 4 , the normalized basis of holomorphic 1-forms Ui can be obtained as follows. Differ- 
entiating (13.17P with respect to a\ 



fi = -i- A ( x d JL. 

J 2m da? J Ai y 



(3.19) 



Comparing with the first equation in fl3.3|) . this means that 



d f dy 



(3.20) 



16 Because of the subtlety mentioned above about how to take 1-cycles that lifts to compact 3-cycles in the 
Calabi-Yau, we should think of the A 1 appearing in (|3.17| e.g. as A 1 = A 1 — A N , where i = 1, . . . , N — 1. For 
simplicity of presentation, we write A 1 as A % . 
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where the total derivative term drj is fixed by requiring that Ui = 0(x 2 )dx as x — > oo. 
Specifically, this leads to drj = d(-xlogy) and Ui is given by 

M , = A ( - logydt) = - "/f^d,. (3.21) 



Although logy may appear problematic because it is not single- valued on the Riemann surface, 
its a 1 derivative is single-valued and does not cause any problem. 

As we discussed in general in section 13.11 turning on noncompact flux breaks M = 2 
supersymmetry to M = 1 by inducing a superpotential. In the present case where the Riemann 
surface is hyperelliptic, we can take {£ m } and {uji} to be the specific ones given in Appendix 
IA.ll As in (13. 8p . the 3-form flux in M.^w reduces to a harmonic 1-form on S sw : 

JV-1 N-l 

m>l i=l i=l 

Under the condition that the compact flux vanishes, (13.111) . this leads to the scalar potential 

(En}. 

We can write the superpotential we are adding to the system in a form that will be useful 
later. By manipulating the quantity Kj n T n appearing in (I3.12p . 

K jn T n = T n I i n = -2T n I x n u = 2T n ^- ( I x n \ogydx 

JBj Joo V a j \Jod 

2T" 3 fl x n +1 dv\ {X2S) 



n + 1 daj V. Joo y 

Here we used (13. 7p . (1A.14[) . and (I3.2ip . By examining (13.121) and (I3.18p . one sees that the 
superpotential is given by: 

W / sw = $^T m Mm+1 , (3.24) 

rn 

where we defined 

' 1 ^m-l \ 1 1 ^m^y 



u m = —— <b x m ~ x \ m = —— (b x m ^. (3.25) 

2mm Jcc 2mm J^ y 

So far everything was about geometry. Now let us turn to the gauge theory interpretation 
of these. As we mentioned above, the local Calabi-Yau geometry (13.131) without flux realizes 
Af = 2 Seiberg-Witten theory, with the hyperelliptic curve (I3.14J) identified with the J\f = 2 
curve of gauge theory. The special coordinates a 1 defined in (13.171) correspond to the U(l) 
adjoint scalars in the IR and parametrize the Coulomb moduli space. The superpotential 
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(j3.24p also has a simple gauge theory interpretation. To see it, we need the relation between 
the vev of the adjoint scalar $ and the curve £sw> given by [54,55]: 

dx \ dy P'n(x) 

tr > = — = NK ' ^=dv. 3.26 

x-$/ y ^P N {xf-Ah? N 

In other words, u m defined geometrically in (13.251) has an interpretation in gauge theory as 
follows: 

« m = -(tr$ m ). (3.27) 

m 

From this, one immediately sees that the superpotential (j3.24p can be written as 

w ™ = E ^qri tr $m+1 = tr [^($)L ( 3 - 28 ) 

m 

where we defined 

™ 

w(x) = J2 — xm+1 - ( 3 - 29 ) 

m 

In (13.281) . $ is understood as the chiral superfield whose lowest component is the adjoint 
scalar. 

Therefore, the J\f = 2 gauge theory perturbed by the single-trace superpotential (13.281) 
corresponds to the geometry ( 13. 13f) with the flux g obeying the following asymptotic boundary 
condition: 

g ~ J2 ™T m x m ~ l dx = W"{x)dx, (3.30) 

in 

where we used (lA.lip . Note that this equivalence holds for any configurations, supersym- 
metric or nonsupersymmetric, because we have shown the equality of the full off-shell scalar 
potential. The perturbed M = 2 theory is precisely the system which was shown in [24, 371 
to have nonsupersymmetric metastable vacua if the superpotential is chosen appropriately^ 
Therefore, it tautologically follows that the IIB Seiberg-Witten geometry with flux at infinity 
also has metastable vacua, if we tune the parameters T m appropriately. 

As we mentioned above, this IIB Seiberg-Witten geometry is dual to a IIA brane con- 
figuration of NS5-branes and D4-branes which can be lifted to an M5-brane configuration. 
In [56], it was shown that superpotential perturbation corresponds in the M-theory setup to 



17 It was shown in [24] to be possible to create metastable vacua by a single-trace superpotential of the form 
(I3.28[) at any point in the Coulomb moduli space for SU(2) and at least at the origin of the moduli space for 

SU(N). 
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"curving" the M = 2 configuration of the M5-brane at infinity in a way specified by the super- 
potential. The metastable gauge theory configuration of [24,37] was realized as a metastable 
M5-brane configuration and its local stability was given a geometrical interpretation. The 
above proof of (13.281) is exactly in parallel to the one given in [56] for the M-theory system. In 
passing, it is also worth mentioning that the M-theory analysis of [56] revealed that at strong 
coupling the nonsupersymmetric configuration "backreacts" on the boundary condition and 
it is no longer consistent to impose a holomorphic boundary condition specified by a holo- 
morphic superpotential, which is in accord with [17]. Therefore, also in the IIB flux setting, 
it is expected that if we go beyond the approximation that the flux does not backreact on 
the background metric, nonsupersymmetric flux configurations will backreact and it will be 
impossible to impose a holomorphic boundary condition of the type ( 13. 22ft . 

Although we do not do it in the present paper, from the viewpoint of flux compactification, 
it is a natural generalization to consider fluxes through the compact 3-cycles. Such flux will 
make additional contribution to the superpotential of the form e^a 1 + m l Fi (see 12.5. IT) . On 
the gauge theory side, in the Seiberg-Witten theory, this can be interpreted as perturbation 
one adds at the far IR, but its UV interpretation is not clear [41]. 

3.1.2 Dijkgraaf-Vafa (CIV-DV) Geometries 

Another example of geometries of the type (13. ip is type IIB on 

M m : uv -F DV (x,y) = 0, x,yeC, (3.31) 
where the underlying Riemann surface Sdv is a hyperelliptic curve 

Sdv : F DY (x,y) = w 2 - [P n (x) 2 - / n _i(s)] = (3.32) 

and P n (x) and f n -i(x) are polynomials of degree n and n — 1, respectively. If we write 

n-l 

f n . 1 {x) = Y,b i x i , (3.33) 

then the coefficients of P n (x) as well as 6 n _i are nonnormalizable and fixed 18 !, while 6j, i = 
0, . . . , n — 2 are normalizable complex structure moduli. The holomorphic 3-form is Qdv = 
— A dx A dy which reduces to 

Adv = xdy (3.34) 



18 More precisely, 6„_i is log-normalizable and can be treated as a variable modulus if one wishes, but in 
the present paper we treat it as a non-dynamical parameter. 
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on the Riemann surface S DV . The geometry (j3.3ip was studied by Cachazo, Intriligator 
and Vafa (CIV) [57] (see also [58]) in the context of large N transition [59,60] and further 
generalized in [61,62]. The Dijkgraaf-Vafa (DV) conjecture [63-65] was also based on the 
same geometry. We will refer to this geometry as the CIV-DV geometry (13.311) or as the 
Dijkgraaf-Vafa geometry henceforth. 

The structure of the underlying hyperelliptic Riemann surface Sdv (13.321) is similar to the 
Seiberg-Witten case (13.141) ; S DV is a genus n — 1 surface with two punctures at infinity. If we 
represent Edv as a two-sheeted x-plane branched over 2n points, those infinities correspond to 
x = oo on the two sheets. The coefficients of P n (x), which are nonnormalizable, determine the 
position of the n cuts on the x-plane, while the coefficients of f n -i{%), which are normalizable, 
are related to the sizes of the cuts. The first homology Hi(Sdv) is spanned by n — 1 pairs 
of compact A- and -B-cycles (A 4 , Bj), i,j = 1, . . . ,n — 1 with in addition a closed cycle A°° 
around one of the infinities which is dual to the noncompact -B-cycle B^ connecting two 
infinities. Because x,y G C, compact A- and -B-cycles on S^v are all contractible in the 
x,y-plane and hence all compact 1-cycles on S DV lifts to 3-cycles in A-^dv with S 3 topology. 

The special coordinates (J3.4J1 in this case is conventionally denoted by S l , Uf. 

S l = - — : / A D v, n, = - — : / A D v, (3.35) 

2m J Ai 2m J B , 

for which, as in (I3.6p . the moduli space metric takes the special form: 

d8 2 = (lmr i j)dS i d&. (3.36) 

Just as in (13.211) . we can write the basis of holomorphic 1-forms iOi using S l as: 

Ui = 7)Q~i(~ y ^ = o /oT^ 1 , \ dX - ( 3 - 3? ) 

dS 2^P n {x) 2 - f n -l(x) 

Adding flux at infinity and breaking Af = 2 supersymmetry to M = 1 go just as in the 
Seiberg-Witten case. The Riemann surface Edv is hyperelliptic and we take {£ m } and {uj,i\ 
to be the ones given in Appendix IA.1I Just like (13. 8p and (I3.22p . the 3-form flux in A^dv 
reduces to a harmonic 1-form on £e>v : 

N-\ N-l 

g = J2 Tm ^ + Yl hiui + Yl ¥w *- ( 3 - 38 ) 

m>l i=l i=l 

Under the condition that the compact flux vanishes (eq. (I3.1ip ). the 1-form (I3.38P leads to 
the scalar potential (13.121) which, just as we derived (13.241) . can be shown to correspond to 
the following superpotential: 

W BV = Y,T m Z m+1 , (3.39) 

in 
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where we defined 

S m = -i- <f x m ^\ DV = -K- <f x m dy. (3.40) 

The 1-form Adv depends on the complex structure moduli S % of the Riemann surface f)3.32|) . 
Therefore, by changing the parameters T m , we can generate a superpotential which is a quite 
general function of S l, s. The OOP mechanism [24] states that, if one tunes superpotential 
appropriately, one can create a metastable vacuum at any point of the M = 2 moduli space. 
Therefore, also for this Dijkgraaf-Vafa geometry, we expect to be able to create metastable 
vacua by appropriately tuning T m , i.e., flux at infinity. Indeed, in the next subsection we will 
demonstrate the existence of metastable vacua in a simple example. 

We have been focusing on the case where there is flux at infinity but there is no flux 
through compact cycles. However, let us digress a little while and think about the case where 
there is flux through compact cycles but there is no flux at infinity. In this case, the IIB system 
has a standard interpretation [57,58,63-65] as describing the IR dynamics of M = 2 SU(N) 
theorjo broken to Af = 1 by a superpotential W = tr[W n ($)], W^{x) = P n (x), with the 
moduli S l identified with glueball fields. More precisely, if there are N l units of flux through 
the cycle A 1 , where N = X^iV 4 , then the system corresponds to the supersymmetric ground 
state of SU(N) gauge theory broken to \Y[ i SU(N l )~j x U(l) n ~ 1 . It is important to note that 
this equivalence between the Dijkgraaf-Vafa flux geometry and gauge theory is guaranteed 
to work only for holomorphic dynamics, or for the F-term. On the geometry side, one is 
considering the underlying geometry (j3.3ip determined by P n (x) and small flux perturbation 
on it. On the gauge theory side, this corresponds to the limit of large superpotential, where 
one has no control of the D-term. Therefore, there is no a priori reason to expect that the 
D-term of the Dijkgraaf-Vafa geometry, which governs e.g. existence of nonsupersymmetric 
vacua, and that of gauge theory are the same, even qualitatively. After all, two systems are 
different theories and it is only the holomorphic dynamics that is shared by the two! 20 ! 

Despite such subtlety, it is interesting to ask what is the gauge theory interpretation of 
adding flux at infinity, in addition to flux through compact cycles. It is known that the curve 
(j3.32p is related to the vev in gauge theory as [54, 63-66] : 



W 



2 



327T 2 \ X - $ 



tr — )dx = ydx = \/P n (x) 2 — f n -i(%) dx. (3-41) 



19 This is the case when we treat b n -i as non-dynamical. If we regard this as dynamical, this system realizes 
U(N) theory. 

20 Of course, it is a logical possibility that even the D-terms of the two systems are identical, or closely 
related to each other. 
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where W 2 = W a yV a and W a is the gaugino field. Comparing this with (13.401) . one finds that 
the quantity S m defined geometrically in (I3.40p has the following interpretation: 

' (tr(W 2 $ m_1 )). (3.42) 



"' 32vr 2 
Therefore the superpotential (I3.39J) can be written as 

W ™ = ^E Tmtr [ w2$m ] = 3^2 tr[W 2 M($)], (3.43) 

m 

where we defined 

M(x) = ^T m x m . (3.44) 

m 

Therefore, flux at infinity of the following asymptotic form: 

g ~ J2 mT m x m - 1 dx = M'(x)dx, (3.45) 

in 

corresponds in gauge theory to adding a novel superpotential of the form (13.431) . Again, this 
correspondence must be taken with a grain of salt, since it holds only for holomorphic physics. 
Note also that flux through compact cycles will induce glueball superpotential [57] of the 
form aiS 1 + iVTIj(S') added to ( 13.391) . Because this part does not contain tunable parameters 
such as T m that can be made very small, it is difficult, if not possible, to use the OOP 
mechanism to produce metastable vacua in that case. 

Now let us come back to the main focus of the present paper, the case where there is no 
flux through compact cycles. In this case, we do not have an interpretation of the system 
as such an SU(N) theory described above, simply because N = ^ N l = 0. Below, we take 
the Dijkgraaf-Vafa geometry with flux at infinity and no flux through compact cycles as an 
example, and see that we can generate metastable vacua by adjusting the parameters T m 
using the OOP mechanism outlined in the previous section. 

3.2 Metastable Flux Vacua in CIV-DV Geometries — An Example 

To demonstrate that one can truly realize supersymmetry-breaking via the OOP mechanism 
in type IIB Dijkgraaf-Vafa flux geometries, we turn our attention now to a simple example, 
namely the geometry relevant for SU(2) 

uv - F DV (x, y) = 0, with F DV (x, y) = y 2 - [P 2 (^) 2 - hx - b ] , (3.46) 
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where we choose 

A 2 
P 2 (x) = x 2 -—. (3.47) 

For simplicity, we will impose a Z 2 symmetry on the Calabi-Yau under which x <-► —x, the 
effect of which is to set the log-normalizable modulus b\ to zero 

61 = 0. (3.48) 

As usual, we can focus our attention on the associated Riemann surface, Suv, wmcn in 
this example has genus 1 and is determined by the equation 

F DV (x,y) = y 2 - [P 2 (x) 2 - b ] = 0. (3.49) 

This geometry admits a single dynamical modulus, bo- This, in turn, can locally be traded for 

the special coordinate S 1 which, for notational simplicity, we refer to as S in the remainder 

of this section 

S = S 1 = — I xdy. (3.50) 

2m J A i 

Alternatively, we can parametrize the moduli space by the globally well-defined coordinate 

£ 2 ( 13.401) . which we choose to denote simply by E 

£ = £ 2 = — J> x 2 dy. (3.51) 

4vrz J x=00 

To this geometry, we now consider turning on flux given by 

g = J2T m U + hu + lzJ, (3.52) 

m>l 

where ui is the unique holomorphic 1-form on £dv- As we have seen, this induces a nontrivial 
potential for £ of the form 



- (^) (§) (^) (§) (^> 

where 

^ DV (S 2 ) = ^T m S m+1 (S 2 ). (3.54) 

m 

To engineer a metastable vacuum at a fixed point, YS°\ we need only choose the T m so that 
Wdv(^) is a cubic polynomial in E obtained by truncating the expansion of a Kahler normal 
coordinate associated to S^ ) at cubic order. To determine the requisite T m , we proceed 
in two steps. First, we must determine the relation between S and the higher S m (13.40)) . 
This is rather trivial. Second, however, we must obtain an expression for the Kahler normal 
coordinate associated to a generic point S^ \ This will be slightly messier. 
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3.2.1 Relating E m and E 2 



Evaluating generic E m for Edv is relatively easy to do given the denning equation (13.490 and 
leads to the result 



E^ 
^4q-2 



0. 



E 

?1=0 



T{2q-n + \) /A 
2^{2q-2n + l)\n\ \~2 



2\ 2(q-n)+l 



fc-^Y 

16/ 



2?-l \ b °~Wj 

^2v^F(2g-2n)!n! 




(3.55) 



A 2 
~2~ 



2(q-n) 



*-£)' 



From this, we first see that E is proportional to bo 

E 



bo 
2' 



(3.56) 



More importantly, however, we are also able to immediately read off the degree of each nonzero 
E m when viewed as a polynomial in E 



E4«j ~ 



C(E 



q-ls 



J iq-2 



C(E 



g-i^ 



(3.57) 



Consequently, the lowest m for which E m contains a term proportional to E 9 is m = 4g — 2. 
This means that to introduce terms of order E 3 into V4 / dv(E), it will be necessary to include 
E m up to m = 10, leading to much more singular flux than one might have otherwise thought. 
This is the first example of a general lesson we will have more to say about later, namely that 
when engineering OOP vacua, the requisite noncompactly supported flux can have a large 
degree of divergence which, to the best of our knowledge, is not easy to determine by any 
simple arguments. 

Because one can introduce quadratic (cubic) terms using any E m with 6 < m < 9 (10 < 
m < 13) there is some choice as to which T m we can turn on to achieve a particular desired 
Wdv(^). For the purposes of our example, we will only turn on T 1 , T 5 , and T 9 , thereby 
adding terms proportional to 



E 2 
E 6 

Eio 



s. 

1 

16 

1 



(EA 4 + 8E 2 ) 



(3.58) 



256 



(A 8 E + 48A 4 E 2 + 128E 3 ) . 



29 



3.2.2 Kahler Normal Coordinate for E 

We now proceed to the second step, namely computing the first few terms of the Kahler 
normal coordinate expansion about a generic point £(°) 

z=(p- E<°>) + a 2 (£ - £(°)) 2 + a 3 (£ - E^) 3 + . . . , (3.59) 

where we have implicitly defined the coefficients 



(3.60) 



fl 2 — 2 £E; 

as = Q^dj, (<fe S r| s ) , 
with (yf E ^ the metric associated to the £ coordinate 

£es= ^ Imr, (3.61) 

and r^ s the associated nonvanishing Christoffel symbol. Computations of quantities such 
as dS/dT, in Dijkgraaf-Vafa geometries are often performed using a perturbative expansion 
about the singular point 5* = 0. For engineering OOP type vacua, though, we need to consider 
instead the neighborhood of a generic, nonsingular point So away from S = 0. Fortunately, 
in the simple case of a genus 1 curve, we can actually obtain exact results without too much 
work by taking advantage of the parametric description reviewed in Appendix [B] As described 
there, one finds that both S and E can be expressed directly as functions of r 

^[l 2 rf (3 ' 62) 



A 4 (12p(r/2) 2 -s 2 ) 



(3.63) 



288p(r/2) 2 

where p(z) is the Weierstrass p-function, g 2 is the Weierstrass elliptic invariant appearing in 
the relation 

= ip{zf - g 2 p{z) - &,, (3.64) 



dz 
and ?7i is one of the half-periods of the Weierstrass ^-function 

Vi = C QV (3.65) 

From ( 13.621) and ( 13.631) . we can apply the differentiation formulae of Appendix [B] to write 
both dH/dS and g^ as functions of r 

<9E in A 3 , . 9 ... i 

=» <? S £ = - AV(t/2)Ihit. (3.66) 



95 v^(772) 
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It is now straightforward to determine the coefficients of the Kahler normal coordinate ex- 
pansion (I3.59P in terms of the value of r at S^ ^ 

36p(r/2) 2 (g 2 + 12r^(r/2) - 6p(r/2) 2 ) - 21 ^J 2 ? 



a 2 



A*(g 2 - 3p(r/2) 2 ) (g 2 - I2p(r/2f) 

864p(r/2) 4 



° 3 A* (g 2 - 3p(r/2)2) 2 (g 2 - 15p(r/2) 2 ) 2 
3607rp(r/2) 3 - A8ng 2 p(r/2) 



Imr 

+ hg 2 2 + 9Qr ]l g 2 p(r/2) - Q3g 2 p(r/2) 2 - 720r/ 1 p(r/2) 3 + 252p(r/2) 4 



(3.67) 



3.2.3 Noncompact Flux for Engineering OOP Vacuum 



We are finally ready to explicitly write the noncompact flux needed to engineer an OOP 
vacuum at a generic point E^. In particular, we seek to specify values for the coefficients 
T m which render 

W D v(E) = J2 ^ m S m (E) (3.68) 

in 
equivalent, up to a constant shift, to a truncation of the Kahler normal coordinate expansion 
(13.59P about E^ ^ at order E 3 . Using (13.581) . it is easy to see that the following choice of 
nonzero T m does the job 

r i = 2-f(A* + 16E )+2a 3 g + ^ + 3I-) 

1* = 1|? _ ** ( A 4 + 8So) , (3.69) 



> 



9 ' 

where a 2 and a 3 given by the expressions in (13.67P evaluated at the value of r corresponding 
to E(°). 

These expressions, while nice and exact, are a little cumbersome so let us also consider a 
special case where things simplify. To that end, we try to engineer an OOP vacuum at the 
special point t — i corresponding to a square torus. In this case, several elliptic quantities 
simplify 

ViUi = \ 9s\ T=i = =► 9*\ T=i = Mr/2)\ =i . (3.70) 

The value of E at r = i can be obtained by applying (I3.70p to ([HJ 



A 4 
E„ = -. (3.71) 
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Figure 3: Plot of V(r) in the neighborhood of our engineered OOP minimum at r = i 



This means that the curve ( 13.491) is given at this point by 
2 



y 



A 2 A 4 

4 *-* 2 , LX 

x ~ x + 777 

144 



' ! 'I + 3V2 



A 1 



,r 



1 



4 3^ 



A^ 



(3.72) 



The coefficients ei2 and 03 ( I3.67P appearing in the Kahler Normal Coordinate expansion (13.591) 

simplify to 

9 , , 1080 , „ . 

a>2\ T=i = — and a 3 | T= , ; = -^-5- (3.73) 



^4 Jir=i ^5 

Inserting these into ( 13.59J) . we find that our desired effective superpotential Wdv(^) is given 

by 

siy 12 m«ny" 3 

(3.74) 



rrr ,„s nv, §1E 2 1080S 3 

VyDv(^) = constant + 3L — — h 



A 4 
while plugging into (I3.69J) yields the T m that do the job 



T 2 



885 



t~i(i 



5832 



T 



10 



A 8 
2400 



5A- ' A'" < 375 > 

While metastability of the vacuum at r = i is guaranteed by the OOP procedure, it is also 
gratifying to see it graphically by explicitly plotting the potential near r = i as in figure [3j 

3.3 Degree of Superpotential Required for Metastable Vacua 

As we have seen in the above example, there is an issue about the degree of superpotential 
we have to consider in order to create OOP metastable vacua. In this subsection, we analyze 
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this issue. 

As one can see from (12.341) . (I2.35p . in order to create an OOP vacuum at a specific point 
X(°) in the moduli space, one must be able to adjust the coefficients in the superpotential up 
to cubic terms in AX = X — X(°\ If the dimension of the moduli space is d, this means that 
we generically need to tune 

d(d+l) d(d + l)(d + 2) d ( d+ i)( d + 5 ) 

t/H 1 d= = C d (3.76) 

2 o o 

parameters in the superpotential] 21 ! The last term is subtracting the degrees of freedom to 

choose the vector fcj. 

In the local Calabi-Yau geometries we have been considering, the superpotential is parametrized 

by the coefficients T m . For example, in the Dijkgraaf-Vafa geometry, the superpotential was 

given by (13.391) : 

W DY (S) = J2 T m Z m+ i(S), Z m (S) = — L- I x m dy(S), (3.77) 

m>l J °° 

where we wrote the dependence of S m 's on the moduli S = {^{explicitly. Therefore, if 
E m (5') , s are generic functions of S then, by tuning C n _i parameters 22 ! T 2 , T 3 , . . . T Cn ~ 1+1 , one 
can create a metastable vacuum at a generic point S = S^°\ More precisely, the OOP mech- 
anism requires that, when we expand S m (5)'s around S^ in AS" = S — S^ ' , the coefficients 
of AS, (AS) 2 , (AS) 3 terms are all independent and by taking linear combinations of S m (S')'s 
we can obtain the superpotential (12.351) . 

However, as we saw in the example above, the situation is not generic for small n and 
we need a more detailed analysis about how high degrees one should go, which is done in 
Appendix O The result (eq. ( IC.4J) ) is that, if we would like to make a critical point at a 
generic point in the moduli space, we have to tune on T m at least up to m = m min , wherq 2 

n = 2 : m min = 10, 

n = 3 : mmin = 15, 



S 2 

n° n n 



(3.78) 



n > 4 : m min = — + — + - + 3. 
o 2 3 



There is certain genericily assumption on the dependence of S m on the moduli (see Appendix 
ICl) . and hence the actual degree m one must consider can be larger than the one given above. 



21 For having a metastable vacuum, the superpotential does not have to be exactly the same as the ones 
given in (|2.35[) : if the coefficients are very close to the ones given in (|2.34|) . (|2.35|) . one still expect to have 
metastable vacua. However, this does not generically affect the number of parameters we need to tune. 

22 Note that the number of moduli is n — 1 because we are treating 6 n _i dynamical. 

23 This result is for the case where 6 n _i is treated nondynamical. For the result in the case where b n -\ is 
regarded as a modulus, see Appendix ICl 
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Therefore, in order to stabilize metastable vacua made of n cuts by the OOP superpo- 
tential, we have to consider E m 's up to rather high degree m min given by (I3.78J) at least. 
Because the degree m corresponds to the order of divergence of the flux at infinity (£ m ) , the 
noncompactly supported flux must diverge at infinity at the corresponding speed. 

4 Factorization 

4.1 The Basic Idea 

In the previous sections we described how we can generate a supersymmetry breaking potential 
for the complex structure moduli of a local Calabi-Yau singularity by the introduction of 3- 
form flux which has support at infinity. Allowing flux with noncompact support may lead 
to various conceptual difficulties, such as the divergence of the total energy density. To 
clarify these difficulties we would like to sketch how such a system can be interpreted as an 
approximation of a larger Calabi-Yau threefold with flux of compact support in a certain 
factorization limit. 

More precisely, we start with a Calabi-Yau with a subset of cycles pierced by usual 3-form 
flux of compact support. In another region of the manifold we have a second subset of cycles. 
The flux from the first cycles will generate a potential for the complex structure moduli of 
the second set. In the limit where the cycles are separated by a large distance, and where we 
zoom in towards the second set, the flux from the first subset will look as if it is coming from 
"infinity"! 24 !. In this sense, the noncompact setup considered in the previous section can be 
considered as a small part of a larger Calabi-Yau with compactly supported flux. 

In this section we would like to understand this embedding into a bigger Calabi-Yau in 
more detail. Our goal is to see how the potential (12.301) arises starting from the standard 
Gukov-Vafa-Witten superpotential for 3-form flux in the larger Calabi-Yau. 

For simplicity we will work with a noncompact Calabi-Yau A4, 

M: uv-F(x,y) = 0, (4.1) 

which is based on a Riemann surface £ given by F(x, y) = 0. As we explained before the 
complex parameters entering the defining equation of the Riemann surface correspond to 
complex structure moduli of the Calabi-Yau. Some of them are non-normalizable and can be 
considered as external parameters. We want to tune these parameters to approach the limit 
where the surface S factorizes into two surfaces E^ and £# connected by long tubes. This 
factorization lifts to the entire Calabi-Yau M. and divides it into two regions M.l and A4r 



As we will see in more detail later, we also have to scale the flux in an appropriate way. 
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that are widely separated. We introduce 3-form flux G3 of compact support on the 3-cycles 
of A4r. The superpotential and scalar potential are given by 

W= fG 3 Att and V = G a d I WdJW, (4.2) 

where the indices /, J run over all complex structure moduli of the total threefold M.. Using 
the properties of the Kahler metric Gj-j in the factorization limit we show that the part of the 
potential (j4.2p which depends on the complex structure moduli of A4l is of the form f)2.30|) . 
Furthermore, we find an understanding of the effective value of the parameters T m . 

4.2 Geometry of Factorization 

In this subsection we study the degeneration of a Riemann surface £ into two components 
Tjl and £#, depicted in figure l4l 25 l In this factorization data of the full Riemann surface is 
expressed in terms of the complex structure of the individual surfaces. It is well known that 
in the limit where the length of the tubes L = 1/e goes to infinity the period matrix of the 
full surface becomes block diagonal 

- LL \ 
t rr)+0 (e) . (4.3) 

While the off-diagonal components t lr go to zero in the factorization limit, their sub leading 
behavior is quite important in our analysis since it expresses the weak interaction between 
the two sectors. The period matrix t lr can be computed systematically in an expansion in 
e from data on each of the two surfaces as we explain below. 






Figure 4: Two conformally equivalent ways of viewing the factorization of a Rie- 
mann surface into two parts. Physically though, we should distinguish both points 
of view, since particle masses depend on the size of the cycles. Because in our 
situation no new massless appears in the factorization limit, the left diagram rep- 
resents our point of view best. 



Technically, we describe the factorization of the Riemann surface with the plumbing fixture 
method [67]. So consider two Riemann surfaces S^ and E# of genus gz and gn respectively. 



25 In general these components could be connected in a non-trivial way. We restrict our computations in 
this section to the case in which they are linked by just one long tube. These should be easily extendible to 
more general cases. 
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n] u a > = ** 2^'-^ T " (4 ' 4) 



On the left surface Ex we have g^ holomorphic differentials Ui, while on the right surface 
E^ similarly g^ holomorphic differentials city. The complex structure of the left surface is 
determined by the periods of the holomorphic differentials 

1 /" ,, 1 

"' 'Bi 

where r^ is the period matrix of E^, and we choose our definitions similarly for the right 
surface. 

The plumbing fixture method works after choosing a puncture P on E^ and P' on E_r. 
It connects the two surfaces by a long tube of length L which is glued onto neighborhoods 
of the punctures P and P' . More precisely, we pick a local holomorphic coordinate z around 
the puncture P such that z(P) = and a holomorphic coordinate z' near P' with z'(P') = 0. 
Then we identify points in these neighborhoods as 

zz' = e. (4.5) 

Now we want to compute the period matrix of the full Riemann surface in terms of complex 
structure data of the two surfaces. For this we need to understand how the differentials u>i and 
ujii extend to well-defined holomorphic differentials on the full surface E = E^ UE^/ ~, where 
~ is the above identification. Let us first consider how to lift the differential u>$. Around the 
puncture P it may be expanded as 



m=l 



J2K? m z m - l dz, (4.6) 



where the functions K[ m are given by ( 1A.7I) . Once we write this in terms of z' we observe 
that, as seen from the right surface, the differential has a Laurent expansion. So Ui will be 
written as a linear combination of the meromorphic differentials £^ of the right surface. A 
meromorphic differential has the following expansion around the puncture 



^-{J^+^n^jdz. (4.7) 

Here we have introduced the functions h^ n , which depend on the complex structure moduli 
of the surface and the position of P. So in general the differential u>i will lift to a differential 
uji on the full surface which can be written as 



OJi 



DC 



^ + > z;<m£m on 



oo m=1 (4.8) 

X] Vimim On E fl . 



v m=l 
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for some coefficients Xi m and yi m . Matching the differential on the two sides we find the 
following conditions 



m N 



x im — 7 .Vvnft'nmi Vvrn ~ I **~im ' / .^w^rvm 1 ■ l^.yj 

n=l \ n=l 

This allows us to compute the cross-period matrix as 

p oo oo n 

t lr - / ,, - S^ k p ' i, - - X^ —k p r- 1 K p> 

'ij' / ^i / j ^i'milim / j ly imy r mn 1 ^j'ni 

Jb,, 1 1 n 



j' m=l m,n=l 

00 € n+l 
^mn == Omn / 'T'^ml "in- 

l=l 



(4-10) 



From this equation we can read off all order e-corrections to the off-diagonal piece of the 
period matrix when a surface E degenerates. 

Also, this procedure gives a clear understanding of the term "flux at infinity". We see that 
the flux at infinity is generated by regular forms on the degenerated surface, and therefore 
will at most have finite order poles at the punctures. 

Notice that for a Calabi-Yau threefold that is based on a Riemann surface, the factorization 
region is described by the deformed conifold geometry 

uv + x 2 + y 2 = e, or equivalently uv + zz' = e. (4-11) 

Usually, this is described as a 3-sphere shrinking to zero-size when e — > 0. However, as for 
the complex 1-dimensional plumbing fixture case we want the two sectors to be far apart 
from each other. Therefore we consider the conformally equivalent setup where the 3-sphere 
is scaled to be of finite size, while the transverse directions are made very large. The finite 
size three-sphere reduces to the cross-section of the tube on the left in figure HI whereas the 
transverse directions reduce to the tube-length. 



To describe the left and right neighborhoods of the degeneration, we can fix x = \Je — y 2 — uv 
on the left and x = — ye — y 2 — uv on the right. In the limit that e — > these neighborhoods 
will not just intersect in a point, but in the divisor uv + y 2 = 0. This is the region where 
regular forms on the total threefold will develop poles when the degeneration starts. 

4.3 Dynamics 

Now we consider turning on flux on the threefold. For simplicity we again take a Calabi-Yau 
(14. ID that is based on a factorized Riemann surface. We turn on 3-form flux G% = F 3 — tH 3 
which is only piercing the set of A-cycles corresponding to E#, as can be seen in figure El 
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and write down the corresponding (super) potential. For regularization issues later, we take 
two more punctures on the right surface labeled by ±00 and turn on some flux a through the 
noncompact B^ cycle running from +00 to —00. 





Figure 5: Turning on flux on the right part of the factorized Calabi-Yau. 



A basis of A and B cycles is given by the compact 3-cycles on the left and the right, 
together with the lift .4.°° of the A-cycle enclosing +00 and B^. So the flux is determined by 



G, = 0, 



A' 



G, = W 



A>' 



G, = 0, 



G, = 0, 



8' 



G 3 = 0, 



Gh = a. 



(4.12) 



Let us denote the complex structure moduli and their duals by X 1 and Fj, which are the 
A 1 resp. Bi periods of the holomorphic 3-form Q. Here we use the capital indices / = {i, i', 00} 
to run over both the left and the right sides. Then the GVW superpotential for the complex 
structure moduli is given by 



W 



f G 3 A n = aX°° + J2 N*'Fp, 



(4.13) 



and the corresponding scalar potential by 



V = J2 G LJ d I Wd J W. 
1, j 



(4.14) 



Since X°° corresponds to a log-normalizable period and the derivatives in the above potential 
just correspond to normalizable modes, the a-factor decouples. This shows that 



i,j,k',l' V /LL i,j',k',l' 

At At L.! It \ / RR 



N k 'r^ 



k'i 



Imr 



{W'-rf?) 



LR 



(4.15) 



i',j',k',l' 
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Thus the total potential is the sum of three terms, which we denote in the obvious way by 

V = V 1 + v 2 + v 3 . 

Next we consider what happens in the limit where the distance L between the two sets of 
3-cycles gets very large. As explained before the period matrices t ll and r RR remain of order 
one in this limit and become almost independent of the moduli Xr and Xl, respectively. 

On the other hand, t lr goes to zero which would make the first term V\ in the potential 
vanish in the limit that e — > 0, at least if we don't scale the fluxes N l appropriately. Since 
V\ describes the interaction between the two sides of the Calabi-Yau, we really want to scale 
the fluxes N l to go to infinity in such a way that the term V\ remains finite. 

Then it becomes clear that the term V3 of the potential dominates over the other two 

contribution to V. This implies that in the limit e — ► the term V3 should be minimized first, 

i.e., 

XX7|g? = 0, Vi', (4.16) 

k> 

which is a set of hr equations for the hr moduli x J . The solutions of this system correspond 

to supersymmetric vacua for the 3-cycles on the right side. Once we have fixed all X^ to 

their supersymmetric values X J , we can consider the effect of the backreaction of the right 

side to the left. This is purely expressed through the potential ]/]_, since the term V2 vanishes 

as well at the supersymmetric point. 

So effectively the potential for the complex structure moduli X l L of the left surface is 

V! = £ (Vr£f) (j^Y WW)- ( 4 - 17 ) 

i,j,k',V ^ ' LL 



This may be written as V\ = J2i j ®i Weg(l/ lm.r) l l L djW e s, where we define the effective "su- 
perpotential" for the left complex structure moduli as 

k' 

Comparing with expression (14.10p it is clear that the fluxes on the right should be scaled 
in such a way that the coefficients 

T m = e m J2 Nk ' K 'k> m (4-19) 

remain constant. In that situation the effective superpotential is 

d i Wtf = y %2T m K im (4.20) 

m 

to leading order in e, which is precisely of the form (12.30J) . 
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4.4 Genericity of Potential and Metastable Vacua 

Let us summarize what we have demonstrated so far. We started with a large Calabi-Yau 
that consists of two parts Ml and Mr separated by a large distance, and turned on a large 
3-form flux on one of the sides, say Mr. This flux generates a large potential for the complex 
structure moduli of Mr, which are therefore set to their supersymmetric minima. The flux 
on Mr is also weakly backreacting to the other side Ml, inducing a small superpotential for 
the complex structure moduli of Ml- We computed this superpotential in equations (14.181) 
and (14.201) and found that it is of the form (12.301) . The main point is that the side Ml only 
knows about Mr via the parameters T m given by (I4.19p . 

In this section we discuss two questions. The first to which degree we can tune the 
parameters T m independently. And the second is whether these T m 's can be chosen to realize 
an OOP supersymmetry breaking superpotential. 

As we can see from (14.191) . the values of the parameters T m depend on the fluxes N l on 
the cycles of Mr and also on the value of the (generalized) period matrix K' Vm . The last 
one depends on the choice of the supersymmetric vacuum X J on the right side. For given 
large fluxes N l there is a huge number of supersymmetric vacua, or solutions of (14.161) . with 
different values of A J and consequently of K' Um . The density of such supersymmetric vacua 
over the complex structure moduli space of Mr has been studied before [6,68-71], and it is 
believed that the vacua become dense in the moduli space in the limit where the fluxes are 
very large. 

The coefficients K' Vm are holomorphic functions over the complex structure moduli space 
of Mr. So naively one would conclude that when the dimension of this moduli space is large 
enough, meaning that the number of 3-cycles in Mr is large, we can always find supersym- 
metric points where the i^'/ m 's have the desired values. However the functions K' Vm are not 
"generic" and there may be relations between them which affect the naive counting. We have 
not analyzed this problem in detail but we think the following statement is true. Any number 
of the T m 's in the superpotential (14.201) can be tuned by considering a Calabi-Yau whose right 
side Mr has a sufficiently large number of 3-cycles, and there will be some supersymmetric 
vacua with right values of K' Vm to reproduce the desired T m 's to good accuracy. 

This claim is made more intuitive by the following physical interpretation of equation 
(I4.19p . Start by turning on fluxes N l on the cycles of Mr, which is based on the Riemann 



surface E^. When reduced on the Riemann surface the flux looks like the electric field 
produced by a charge in two dimensions. The set of fluxes N l resembles a charge distribution 
on the cycles of the Riemann surface. To compute the field produced by these charges in the 
distant region of the other set of cycles £/,, one has to consider a multipole expansion. Since 
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the matrix K[, m computes the mth multipole expansion of a charge distributed along the Z'th 
cycle, the coefficients T m are exactly the multipole moments of the charge distribution. In 
this formulation our first question reads whether we can arrange a charged distribution to 
have the desired multipole moments given by the coefficients T m . We expect that the answer 
is positive. 

The second question is more subtle. To realize a metastable nonsupersymmetric vacuum 
via the OOP mechanism, one has to tune the superpotential in a way which is determined 
by properties of the Kahler metric at that point. As we saw in section 13721 one has to tune 
the coefficients of the effective superpotential only up to cubic order in an expansion around 
the candidate metastable point. Since we have a very large number of parameters T m at 
our disposal it seems that generically we should be able to tune them to generate metastable 
vacua at most points on the moduli space. However we do not have a proof of this statement 
and it is possible that various relations between the period matrices and the Kahler metric 
invalidate the naive counting 2 6 I. 

5 Factorization II: An Example 

In the previous section, we argued, based on the factorization of the Riemann surface and 
Calabi-Yau, that it is possible to embed the nonsupersymmetric metastable vacua we found 
in [3] in a "larger" Calabi-Yau, the idea being that the flux threading compact cycles on one 
side of the Calabi-Yau looks like flux coming from infinity from the viewpoint of the other side 
of the Calabi-Yau. In this section, we will discuss the Dijkgraaf-Vafa geometries of subsection 

EX2 

n 

S D V : V 2 = Pn(x) 2 - fn-l(x), P n (x) = JJ(x - «/), (5.1) 

7=1 

as an example where our proposal can in principle be implemented, and make some steps 
towards actually confirming our proposal. 

5.1 Factorization Limit in Practice 

As explained in 13.1.21 aj are non-normalizable parameters which represent the positions 
of the cuts on the x-plane, while the coefficients in f n -i(x), or equivalently variables S 1 



26 This question is similar to whether one can realize the OOP mechanism with a single trace superpotential 
for the adjoint scalar in an SU(N) gauge theory. In [24] it was demonstrated that for SU(2) a metastable 
vacuum can be generated anywhere on the moduli space by a single trace superpotential, and for SU(N) at 
the center of the moduli space. It was not fully analyzed whether this is possible in generality. 
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defined in ( 13.351) . are normalizable (or at least log-normalizable) and hence are dynamical 
variables describing the size of those cuts. Therefore, in this Dijkgraaf-Vafa case (j5.1| . otj 
are the parameters we want to adjust in order to approach the factorization limit where Sqv 
degenerates into two subsectors. 

So, what we should do is clear: we divide the n cuts into two parts as n = til + Ur, the 
ones on the left indexed by i and on the right by i', and send these two groups apart from 
each other by a large factor L = 1/e so that 

<Xi-a v = 0(L) (whenL-^oo). (5.2) 

In the L — > oo limit, the left and right sides will be very far apart and the factorization we 
discussed in the previous section must be achieved. For example, the period matrix of the 
total Riemann surface must diagonalize as in (14.31) up to 1/L correction. 

There is one thing we should be careful about when taking the L — » oo limit. If we try 
to separate the two sets of cuts by naively taking the typical difference between on and a^ to 
be of order L while keeping the size of the cuts fixed, then a simple estimate of the scaling 
of Si , Sft using (I3.35P shows that the physical size of the 3-cycles in the Calabi-Yau blows 
up. What we want instead is to end up with two sets of 3-cycles of finite size, separated by 
a large distance, so that we are left with nontrivial dynamics of S^, Sfi. To achieve this we 
must also scale the size of the cuts, as we send L —> oo. Let xl and xr be local coordinates 
in the left and right sectors, respectively, and set 

x L = L r x L , x R = L r 'x R , (5.3) 

where 

' HL (5.4) 



n L + 1 n R + 1 

Then, from (13.351) . it is not difficult to see that we can keep S^, Sft finite if we keep xl, xr 
finite while taking the L —>■ oo limit. A similar rescaling of local coordinates must be also 
necessary when taking a factorization limit in any other examples than (15. ip . 

5.2 Computation of Period Matrix 

In the Dijkgraaf-Vafa geometry (15. ip . the period matrix is given by 

Here, JF is the B-model prepotential, which by the Dijkgraaf-Vafa relation [63, 72] is related 
to matrix models. The precise way to scale various quantities to take the factorization limit 
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being understood from subsection 15. 1[ it is in principle possible to confirm our proposal for the 
Dijkgraaf-Vafa geometry using (15.51) . For doing that, it is important to be able to compute the 
prepotential JF for a large number of cuts n. The results from section [3] show that generating 
a metastable vacuum requires quite a lot of coefficients T m . Since we roughly need the same 
number of cuts on the right as the number of tuned S m 's on the left, the total Riemann 
surface must have quite a large number of cuts. So, in this subsection we will explain the way 
to compute Tq and thus tjj for an arbitrary n. 

For Dijkgraaf-Vafa geometries (15 .ip the prepotential JF may in fact be computed for any 
number of cuts n in a number of ways. The most direct way is evaluating the period integrals 
on the hyperelliptic curve. This has been done up to cubic order in S 1 in [73]. Duality with 
a U(N) matrix model [63,72] 



cxp 



.9=0 



d N $ exp 



trW{<$>) 



(5.6) 



where the matrix model action is given by 

W'ix) = PJx) 



i=i 



X 



Oil) 



(5.7) 



makes this computation quite a bit simpler. Let us quickly show this argument [72]. 

The field $ is an iV x N matrix. Say N 1 eigenvalues of $ are placed at the critical point 
x = ai and divide the matrix $ into N 1 x N J blocks $/j, where Y^i=\ N 1 — N. One can go 
to the gauge $/j = for I ^ J by introducing fermionic ghosts in the matrix model action. 
This produces the following extra term in the action, where $/ = $//: 



W ghost = ^tiiBj&Cu + Cj&rBu). 

1+3 



(5.1 



To write down Feynman diagrams, we expand $/ around x = oej as $/ = aj + 0/. A 
Taylor series of W(Q>i) = W(ati + 4>i) around ai yields the propagator and p-vertices for (f>j. 
In particular, this shows that the propagator for 0/ is given by 

1 1 



iiVi 



) 



W"{ai) A/ (5 ' 9) 

where A/ = W"(ai) = YYj^i a i3- Moreover, expanding the ghost action determines the ghost 
propagator to be 

1 



(BjrCjj) 



OLIJ 



(5.10) 
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j* SfSj 




+ rn j- 



4 SjSjS k 



Figure 6: The contribution to J^o,3 given in terms of matrix diagrams. Gray double 
lines represent 0/ fields, while black-and-gray double lines represent BC ghosts. 



and gives the Yukawa interactions between <ftj, Bjj and Cij. 

The contribution to the prepotential JF of order three in the S I: s is given by planar 
diagrams with three holes. Writing down the expressions gi$ and 5^4 in terms of en's and A's 
shows that 

n n n 

^0,3 = ^J UiSf + ^ Ul;jSjSj + ^J UijkSiSjSk, 

1+3 



(5.11) 



where 



1=1 



UI= 3 



KJ<K 



y, 1 1 y^ 1_ 



afrAj 4A 7 



Ul;J 



3 2 

+ 



afjA/ afjAj ajjAj r ^f , a IK 



£^ 



and 



K+I,J 



UlJK = 4 



1 1 

+ ^ + 



1 



a IJ a IK A I ajja JK Aj a KI a KJ A K 



In appendix 0, we discuss the generalization of this result to higher order in S 1 . In 
particular, we compute JF up to S 5 terms. 



5.3 Scaling of Period Matrix 



The method explained in subsection 15.21 allows one in principle to compute the period matrix 
to any order in S 1 for general Dijkgraaf-Vafa curves (j5.ip . Then the factorization limit can 
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be achieved simply by taking the L — > oo limit of the result and one can start looking for 
metastable vacua. In this subsection, as a step towards it, let us pursue a more modest goal 
of seeing the factorized behavior of the period matrix, (14. 3p . 

The form of the scaling can be elegantly derived for any possible contributing matrix model 
diagram to JF . First note that A« scales as L 2r as L — > oo, and A# as L 2r . All propagators 
with indices from either side of the surface have an expansion in terms of otu's and A/'s, and 
thus a scaling in L which is easy to determine. The total scaling of a planar diagram with an 
arbitrary number of these elements turns out to depend just on the number of ghost vertices 
that connect the left side to the right side. It is given by 

where N^i is the number of ghost vertices with external ghost lines indexed by (i,i') and the 
external 0-line by (i,i). Note that in deriving this we assumed the scaling (15. 3p and thus 
S l L , S l R are of order one. 

This shows that a diagram with only indices on the left (or on the right) will be of order 1 
in L. Since such diagrams contribute to the period matrix r^ (or t^y), so this shows that the 
period matrix is of order 1 in L, with corrections in 1/L from diagrams that contain at least 
two loops indexed by i and j. On the other hand, the off-diagonal pieces of the period matrix 
Tu> and Ti'i contain at least one ghost cross- vertex with indices % and i' . These parts will 
therefore scale at least as 1/L. In particular, for large L the properties of the full Riemann 
surface S are determined by those of the two factors £l, T,r, and the period matrix tij indeed 
diagonalizes as in (14.31) . 

Having checked the diagonalization (14.31) . the problem of actually finding an example of a 
metastable vacuum then just amounts to solving equation (I4.19p together with (I4.16P using 
the data from matrix model, for T m giving a metastable vacuum. Solving these equations is 
nontrivial, since the relation between the flux parameters N % on the right and the coefficients 
in the superpotential T m we want on the left are non-linear, although we expect that the 
solutions do exist by the multipole argument we gave in section HI We leave matrix model 
computations up to requisite orders as well as finding the actual metastable vacua by solving 
those equations for the future work. 

6 Conclusion and Generalizations 

Summarizing, we found that turning on flux with support at infinity in local Calabi-Yau in 
type IIB induces superpotential for the moduli in the local Calabi-Yau, thus breaking M = 2 

45 



of the Calabi-Yau compactification down to M = 2. Then we demonstrated that one can 
create metastable vacua by tuning the flux at infinity using the OOP mechanism, using a 
Dijkgraaf-Vafa (CIV-DV) geometry as a primary example. The metastable vacua known to 
exist [24,37] in perturbed Seiberg-Witten theory can also be understood in terms of metastable 
flux configuration. 

Flux diverging at infinity may appear problematic, but in reality a local Calabi-Yau must 
be regarded as a local approximation of a larger compact Calabi-Yau and the flux at infinity 
has a natural interpretation there; there is flux floating around in the rest of the Calabi-Yau, 
which "leaks" into our local Calabi-Yau and just appear to be coming in from infinity. This, 
furthermore, motivates a more natural setting to realize metastable flux vacua: in a part, say 
on the right side, of the full Calabi-Yau Ai, there are some 3-cycles threaded by flux (and 
possibly O-planes to cancel net charge if M. is compact) and on the left side there are some 3- 
cycles without flux through them. If the distance between the left and right sectors is large, the 
full Calabi-Yau M. factorizes into an almost decoupled system of A4l and Mr, and the flux 
in Mr appears to be flux at infinity from the viewpoint of M. L and induces superpotential in 
M.L- By adjusting the number of fluxes in Mr, we can tune the superpotential and generate 
metastable vacua in Ml- This is a very well controlled setting to analyze flux vacua, which 
may shed light on the structure of the nonsupersymmetric landscape of string vacua. We 
also made some steps toward actually embedding metastable vacua in a larger Calabi-Yau as 
sketched above in the case of Dijkgraaf-Vafa geometry by computing certain matrix model 
amplitudes. Actually finding explicit vacua along that line is an interesting open problem. 

Note that we needed just two main ingredients to achieve this result. The OOP mecha- 
nism requires that the complex structure moduli space is special Kahler, and it is important 
that a superpotential for flux is very much controllable by tuning the flux, such as the Gukov- 
Vafa-Witten superpotential. This means that we can generalize the above story to any setting 
which fulfills these two requirements. Other possibilities therefore include M-theory and F- 
theory on Calabi-Yau fourfolds [28,74]. Let us finish by saying a few words on these two setups. 

Compactifying M-theory on a Calabi-Yau fourfold Ma with fluxes yields a three-dimensional 
low energy theory with 4 supercharges. The complex structure moduli of the Calabi-Yau are 
part of the chiral supermultiplets and are described by variations of the holomorphic (4, 0)- 
form Q. In the local limit where the fourfold becomes noncompact, the Kahler potential on 
the moduli space is given by 



K= / fiAfi, (6.1) 

so that the metric on the moduli space is indeed special Kahler. Moreover, it is well-known 
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that the complex moduli may be stabilized by turning on 4- form flux F 4 , which introduces 
the superpotential 

W= f F 4 AQ. (6.2) 

JMa 

The condition for unbroken supersymmetry is W = dW = 0, so that F± has to be a (2, 2)-form. 
Stabilizing the Kahler moduli as well requires that the flux is primitive under the Lefschetz 
decomposition (and in particular self-dual). Turning on primitive (2, 2) flux on some compact 
4-cycles, we can now follow an equivalent procedure as in IIB. 

M-theory compactified on .M 4 is equivalent to compactifying F-theory on M. A x S 1 , at 
least if .M 4 is an elliptically fibered Calabi-Yau. This leads to a four-dimensional space-time 
with 4 supercharges. So again, the Kahler potential is given by (16.11) . and the flux F 4 is a 
primitive (2, 2)-form. The relation with IIB consistently reduces _F 4 to a harmonic (2, l)-flux 
G3. The extra seven-branes that must be inserted in IIB when reducing over a singular T 2 
do not contribute to the superpotential and thus don't play an important role here. 

In particular, consider as an example the local Calabi-Yau fourfold 

u 2 + v 2 + w 2 + F{x,y) = 0, (6.3) 

where all variables are C (or C*) valued, and F(x, y) defines a smooth curve in the x, y-pl&ne. 
Its holomorphic four-form is given by 

_ du A dv ,„ .. 

Q = AdxA dy. (6.4) 

w 

The u, v, w-fiber defines a two-sphere over each point in the x, y-plane, which shrinks to zero- 
size over the curve F(x, y) = Q 27 l . Four-cycles can be constructed as an S 2 fibration over some 
disk D ending on the curve and have the topology of a four-sphere (when x and y E C). 
Notice that the intersection lattice is symmetric now and not simply symplectic anymore, so 
that the bilinear identity takes a more complicated form. However, like in the threefold case 
all relevant quantities reduce to the Riemann surface, and the analysis is similar as before. 
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A Some Basic Results on Riemann Surfaces 

In this appendix we summarize some basic properties of Riemann surfaces [75] . 

A compact Riemann surface S 9 is a one-dimensional compact complex manifold and its 
topology is completely characterized by its genus g. The middle cohomology group has 
dim H 1 (T lg ) = 1g. The intersection form on Hi(E g ,1i) is antisymmetric and by Poincare 
duality unimodular, meaning that we can pick a basis of 1-cycles A 1 , Bj with intersection: 

A i n A j = 0, B % n Bj = 0, A i n Bj = 8), i,j = l,...,g. (A.l) 

Such a basis is unique up to a symplectic transformation in Sp(2g, Z). 

S g has a complex structure moduli space M. g with dim.M 5 = 3g — 3, g > 2. 

A 1-form wona Riemann surface is called a holomorphic differential if in a local coordinate 
patch it has the form: 

uo = f(z)dz, f(z) : holomorphic. (A. 2) 

We will also consider meromorphic differentials, for which we allow the function f(z) to have 
poles at certain points on the surface. Now we present a standard basis for holomorphic and 
meromorphic differentials on a general Riemann surface: 

Holomorphic differentials^] uji: Once we pick a symplectic basis of one-cycles, there is a 
canonical basis of holomorphic differentials u>j, i — 1, ..,g, with the following periods: 

1 I xi l 

x ujj = 0,, — 

'A* 27T* J B 



2mj« Ui = * M < r- U > = T * (A ' 3) 



The (symmetric) matrix t%j is the period matrix of the surface, which depends on the complex 



structure of S fl . 



Meromorphic differentials of the second kindj^l £,^>i : These are characterized by a 
point P on the surface where the differential has a pole of order m + 1 with m > 1. They 



28 These are also called meromorphic differentials of the first kind. 

29 A more common notation in the literature for meromorphic differentials of the second and third kinds is 



gK7^ and d£l ' . 
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are normalized so that in local complex coordinates z where z(P) = they have the Laurent 
expansion: 



dz 

£m ~ m ~^i + regular. (A.4) 



Meromorphic differentials of the third kind, £ ' : characterized by two points P, P', 
where the differential has first order poles with opposite residues. Around P we have: 

£ ' ~ h regular (A. 5) 

and similarly around P' with the opposite sign. 

Notice that we can always shift a meromorphic differential by a holomorphic differential 
without changing the singular part of the Laurent expansions (1A.4I) . (1A.5I) . We can eliminate 
this ambiguity by demanding that the A periods of the meromorphic differentials vanish: 

C = 0. (A.6) 

A i 

In general, it is not possible to simultaneously set the B periods to zero. Instead we have: 



where the matrix Kf m depends on the complex structure moduli of the Riemann surface and 
the position of the puncture P. 

A.l Hyperelliptic Case 

Let us consider the case where E ff is hyperelliptic. For example, the curve appearing in the 
Dijkgraaf-Vafa case, (I3.32p . can be written as: 

n 

y 2 = P n (x) 2 - /„_i(x), P n (x) = l[(x - a^. (A.8) 

%=x 

This curve can be regarded as a two-sheeted x-plane with n cuts and two punctures, the latter 
corresponding to infinities on the two x-planes. Let us denote these points by oo and oo. 
A basis of holomorphic differentials Ui, i = 1, . . . ,n — 1 can be constructed by 

= 9M dx = 3M =dx, (A.9) 

V V P n(*) 2 -fn-l(x) 



where Qi(x) is a polynomial of degree up to n — 2 chosen so that (1A.3J) holds. Note that this 
u>i goes as ~ 0{x~ 2 )dx as x — > cxd,oo, which means that this is regular at x — oo,oo. 
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In the hyperelliptic case, it is convenient to take the meromorphic differentials of the 
second kind, £ m , as 



6. = ^<fc 



^m\p") 



y ^Pn{x) 2 - /„_il 



zdx, 



m > 1. 



x) 



(A.10) 



Here, R m {x) = mx m+n 1 + . . . is a polynomial and the coefficients of x m+n : 
chosen so that 



£ m = ± [mx m 1 + 0(x 2 )] dx, x ~ oo, oo 



(A.n; 



is satisfied. This condition is similar to (1A.4I) . but this £ m has poles at two points, x = 00,00, 



instead of one. The coefficients of x n 2 , . . . , x° are chosen so that (IA.6I) is satisfied 



The meromorphic differentials of the third kind, £0, can be defined likewise using a poly- 
nomial Rn(x) = x n ~ 1 + . . . , where the coefficients are chosen so that 



£ = ^dx = ± 

y 



- + 0{x- 2 ) 

X 



dx, x ~ 00, 00 



(A.12) 



holds and (1A.6J) is satisfied. 

Let us derive a formula that will be useful in the main text. By expanding the right hand 
side of the trivial identity = J s u>i A £ m by the Riemann bilinear identity, one finds 



= ^ / U{ \ ira - \ im\ Wi I + ^2 fUJid 1 £ m 
i V AJ jB i jAJ jB i J p=oc,So^ 

= K im + J2 fuid-^m- (A. 13) 

p=oo,oo " 

Because the behaviors of Ui, £ m at x = 00 is the same as those at x = 00 up to a sign, 



K, 



p=oo,oo 



^ f ^i d l £ m = -2 (b LUid x £. 



-2 <b x m u % 
'00 



(A.14) 



B Parametric Representation of Genus 1 Curves and 
Sample Computations 

In this appendix, we review the parametric representation of the genus 1 Riemann surface 
Sdv of section 13.21 defined by 



= F BY (x,y)=y 2 - 



2 A 1 / 

x - — ) -b 



(B.l) 
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(a) Depiction of Sdv a s double cover of the 
x-plane with compact A and B cycles indi- 
cated 



(b) Depiction of Edv &s fundamental do- 
main in z-plane with the corresponding 
A and B cycles indicated 



and its application to obtaining some of the results used therein. In particular, we think 
of Spjv as a copy of the standard fundamental domain with two marked points, a\ and a^ 



corresponding to the points at infinity on the two sheets. In figures 7(a) and 7(b) , we depict 
both the standard visualization of £dv as a double-sheeted cover of the x-plane as well as the 
parametric one, identifying the standard A and B cycles in the former and their realization 
in the latter. 



B.l Building Blocks 

The embedding of S DV into xy space is obtained by specifying functions x(z) and y(z) which 
satisfy (IB.lj) . The basic building blocks that we use to construct x(z) and y(z) are Janik's 
functions Fi(z) [76] 

4- (B.2) 



and their derivatives 



Fi(z) = \n6(z — di — t) t 

9 N " 



F^(z) 



dz 



FAz) 



(B.3) 



A detailed description of these functions, their properties, and several sample computations 
can be found in Appendix C of [77]. For now, we simply note a few elementary facts. First, 
we point out that F{(z) introduces a branch point at a, while F^ introduces a pole of order n. 
For n > 2 these functions are elliptic while for n = 0, 1 they have the following monodromies 



Fi(z + 1) 

Fi{z + t) 
if^ + 1) 
F^iz + r) 



Hz) 
Fi(z) 4 

F^{z) 
fP(z) 



VK 



27ri(z 



(B.4) 
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It is also useful to record the relation between these functions and the Weierstrass a, (, and 
p functions 

F(z) = \na(z) - %z 2 + mz + ln6»'(f) 
F (1 \z) = CO) + in - 2r] lZ 
F M( z ) = -p( z )-2r ll ( B - 5 ) 



p( n )i 



z =- 



n-1 



where 



d_ 

dz 



* = <( 5 



p(z) n/,2 



Finally, we also recall the differential equation satisfied by p{z) 



dpjz) 

dz 



4p(^) 3 - 92p{z) - g 3 



(B.6) 



(B.7) 



which can also be taken as an implicit definition of the Weierstrass elliptic invariants #2 and 
9s- 

B.2 The Embedding Functions x(z) and y(z) 

Using the building block functions F^ (z) , it is fairly easy to write down embedding functions 
x(z) and y{z) satisfying ( IB. II) . Because x{z) should be locally one-to-one near the marked 
points, we must construct it from functions with single poles, namely F\ and F 2 . On the 
other hand, y(z) ~ x(z) 2 near the marked points so it must contain functions with double 
poles, F{ and F 2 . This leads us to writer 



x{z)=X (Fi 1] - F 2 (1) - [F^{a) - in] 
y(z) = X 2 (if> 



where 



,(2) 



a = Gt2 — «2 



(B.8) 



:b.9) 



Because elliptic functions such as x(z) and y(z) are completely determined by their pole 
structure, it is in fact quite easy to verify that 



Vizf 



x{zf 



A 5 



2\ 2 



bix(z) - b 



(B.10) 



'The constant term that we add to x{z) is added for later convenience. 
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where 

A 2 = 12X 2 p(a) 

&i = -4X 3 p'(a) (B.ll) 

6 = X 4 [I2p(a) 2 - <? 2 ] 

In order to obtain h\ = we set 

a = \ (B.12) 



We can solve for X in (IB. Ill) 



12p(r/2)6 



^ ~ A3 [12p(r/2)s - g 2 ] (R13) 

and use this to eliminate X, thereby obtaining a direct relationship between A, bo, and the 

complex modulus r 

= A'[12 P (r/2)^ fc ) 

144p(r/2) 2 V ; 

B.3 Two Sample Computations 

We now describe two sample computations which illustrate the power of this approach. First, 
we will reproduce a result that is more easily obtained using the explicit representation fIB.lj) . 
Next, we will consider a computation for which the parametric approach is simpler. 
As our first example, let us consider the quantity 

E = ££„/* (B15) 

As we saw in section I3T21 this can be done quite easily using the explicit representation (IB. II) 
with the result 

S = I (B.16) 



We can also write this directly in terms of r using (IB.14J) 



= A*|12 P (T/2)- ga ] 

288p(r/2) 2 v ; 

Let us now see how the result (1B.16[) can be obtained using the parametric representation. 
For this, we write 



Airt 



L £ x{z ?m dz=±fc X* (F« - F« - [**)(«) _ *] ) 2 ( F f) _ F f) * 

(B.18) 
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and expand the integrand near a\. This is straightforward and leads to 



x(z yij^l, "» ( ,. " — -'■'-— + ([z - ai ] )) 

dz l2p{T/2) 2 -g 2 \{z-a 1 f [z - a t f z - a x ) 

(B.19) 
where we have used ( IB. Ill) . The residue appearing in E is now easily read off with the desired 
result 

E = - (B.20) 

Next, let us turn our attention to the computation of 

S= J> ydx (B.21) 

2m J A i 

In the parametric formalism, we write this as 



s _ 1 f ..,_\dx(z) 

~ 2m 



To evaluate this, we will write the integrand as a sum of quasi-elliptic functions and use their 
known monodromies (1B.4J) . Given that the integrand has poles of degree at most 4 with even 



(odd) poles at Oi and a 2 entering with identical (opposite) signs, the general form of this 
expansion is relatively simple 

(if> - Ff ) 2 = a (f« + Ff ) +b (Ff - Ff ) +c (if > + if >) +d (f« - F« + ztt) +e 

(B.23) 
In terms of these expansion coefficients, the monodromies (1B.4I) lead to the simple result 

S = (iird + e) (B.24) 

2m 

In practice, the coefficients a, . . . , e can be by comparing pole structures on the two sides of 
( IB. 231) with the following result when a = r/2 

o = -i 6 = 0, c = 2p(r/2), d = 0, e= ^ + 8^p(r/2) - 4p(r/2) 2 (B.25) 

This means that S is actually given by the relatively simple expression 
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B.4 Some Useful Identities 



Finally, we close this Appendix by listing a few derivative identities that were useful in section 
13.21 First, some derivative identities 

d((z) _ 1 



dr 

dpjz) 

dr 

dr 



2m 

1 

2vri 
1 



-p'(z) + ({z)p(z) 



92Z 

12 



+ 2r h (((z)-zp(z)) 
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2p{zf + C(z)p'(z) -j-^Vi (zp'(z) + 2p(z)) 



2 2 - — 
2m \ Vl 24 



Several of these can be combined in order to derive the additional useful result 



2p(r/2) 2 - f - 4,, iP (t/2) 



We also remind the reader that the partial differential equation 

2 

= Ap{zf - g 2 p{z) 



dp(z) 
dz 
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combined with the fact that 

p'(r/2) = 

implies that the elliptic invariant g^ can be written in terms of g 2 and p(r/2) as 

g z = Ap(r/2f-g 2 p(r/2) 



;b.27) 



(B.28) 

(B.29) 
(B.30) 

(B.31) 



C Independence of E m 's 

In this appendix, we consider the Dijkgraaf-Vafa geometry and examine the dependence of 
S m 's on the moduli S = {S 1 }, or equivalently, on the coefficients b = {6j} of the polynomial 
fn-i{ x ) as defined in (I3.33p . To apply the OOP mechanism and generate a metastable vacuum 
at a point bi = b\ , it is needed that, when we expand E m (6)'s around a point in Abi = bi — b\ , 
the coefficients of Abi, AbiAbj, AbiAbjAbk terms are all independent and by taking linear 
combinations of E m (fe)'s we can obtain the OOP superpotential (12.351) . 

For simplicity, let us first discuss the case where we treat 6 n _i, which is log-normalizable, 
as a dynamical modulus. In this case the number of moduli is n and the number of coefficients 
we would like to tune is, from ( 13.761) . 

n(n + l)(n + 5) 



Cm 



6 



(c.i; 
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Explicitly, E m (6) is given by 



Sm(6) 



1 



lixim 
— Res 

M X=OD 



x m dy = — Res 
m x=oo 



x 



, 2P n {x)p^x) - r n _^) 

2V 'Pn(x) 2 - f n -l(x) 



X 



pi i \ _ Jn-lA-^J 

^ W ~ 2P„(aO 



x 



1 J_ * /n-1 0*0 , 3 ( fn-l{x) 



2 P n (x) 2 8 V P n (x) 2 



5 

+ 777 



/n-i(aj) 



16 V PJx) 2 



+ 



(C.2) 



So, E m (6) are polynomials in Vs. If they are generic polynomials in b with high enough 
degree, then the expansion of S m (fe) around a generic point b^ ' in Ab will have different 
coefficients of Ab, (Ab) 2 , (Ab) 3 terms, for different values of m. If this were the case, then 
the minimum number of S m 's we need to consider would be C n in (1C.1I) . 

However, for small m, E m (p) is not a generic polynomial in bi and we need to be careful. 
From ( 1C.2I) . one can read off the following pattern of dependence of S m 's on fej's: 

• E_ n , . . . , S do not depend on bkS, because the only contributions come from P' n . 

• A term with just one bi (i = 0, . . . n — 1) appears in P! n (f ! n _ 1 / P 2 ) and f n _ x jP n . Such a 
term has degree i — n — 1 in both cases and hence bi first shows up in E n _j. fo n _i appears 
in Ei and 6o appears in E n . 

• The combination 6,6 



-'"./ 



(i,j = 0, . . . n-1) appears in P^f^/P 2 ) 2 and (f n -x/Pn)U'n-x/P%)- 
These terms have degree i + j — 3n — 1 and hence 6j6j first shows up in E3 n _j_j. 6 2 _ x 
appears in E n+2 and 6q appears in E 3n . 

• The combination b^h (i,j, k = 0, . . . n-1) appears in P^(j '' n _jPlf and (f' n _j ' P n )(f' n - 1 l ' Plf ■ 
These terms have degree i+ j + k — 5n — 1 and hence bibjbk first shows up in E5 n _j_j_fc. 

&o appears in E 2n+3 and fejj appears in E 5n . 

From these, we can see that we have to satisfy some requirements. Let us call m of E m 
"order." 

• We need all combinations of Abi Abj Abk, but &jj, which contains (Abo) 3 , does not appear 
until order A n = 5n. 

• The number of possible cubic terms, Abi Abj Ab^, is n(n + l)(n + 2)/6. Cubic terms 
start to appear at order 2n + 3 and therefore, for all possible cubic terms to have 
chance of all showing up in a linear independent way, we need to wait until order 



(2n + 3) + n(n + l)(n + 2)/6 - 1 



B n . 
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• The number of possible quadratic terms and cubic terms is n(n+l)/2+n(n+l)(n+2)/6. 
Quadratic terms start to appear at order n + 2 and cubic terms appear at higher order. 
Therefore, for all possible quadratic and cubic terms to have chance of all showing up 
in a linear independent way, we need to wait until order (n + 2) + n(n + l)/2 + n(n + 
l)(n + 2)/6-l = fl„. 

• From ( 1C.1I) . we need C n independent coefficients. So, we need to wait until at least 
order C n . 

By looking at which of A n , B n , B n , C n is largest for given n, we find that we need S m 's at 
least up to m min , where 



n — 1 > , o7, m i n — o, 

n > o _ m - r _ (n + l)(n + 2)(n + 3) (C3) 



By a similar analysis, if 6 n _i is regarded as a nondynamical parameter, we find the fol- 



lowing: 



= 2 - 


■* "^min = 10, 


= 3 - 


-> ^7imin J-") 


>4 - 


S 2 

n n n 
o 2 3 



(C.4) 



D Prepotential for Dijkgraaf-Vafa (CIV-DV) Geome- 
tries 

In this appendix, we first review different approaches to computing the prepotential T§ for 
the Dijkgraaf-Vafa (CIV-DV) geometries [57,58,63-65] given in eqs. ( 13T3T]) . (I3~32l) : 

uv-F m (x,y) = 0, (D.l) 

F DV (x,y) = w 2 - [P n {xf - /n-i(x)] , (D.2) 

P n {x) = W\x) = g n+ i ]J{x - a { ), (D.3) 

for arbitrary number of cuts n. Moreover, we will present Tq for general n up to S 5 terms. In 
the present paper, the prepotential is used in section [5] to evaluate the period matrix of the 
underlying hyperelliptic Riemann surface. However, the content of this appendix is almost 
independent of the main text and can be read separately. 
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The prepotential is physically important, because by putting fluxes in the Dijkgraaf-Vafa 
geometry one can realize supersymmetric N = 1 U(N) gauge theory, and its glueball super- 
potential which governs low energy dynamics can be computed from the prepotential [57,58]. 
Furthermore, by the Dijkgraaf-Vafa relation [63-65], the prepotential is related to unitary 
matrix models. The relation to matrix models was studied also using supergraphs [78] and 
Konishi anomaly [54]. The same prepotential also underlies the physics of metastable brane- 
antibrane systems studied recently [9,77,79]. 

The first computation of the prepotential was performed in [57] for n = 2 (two cuts) up to 
S 5 terms by directly evaluating period integrals, where S is the glueball. For small values of n, 
the computation of JF up to several orders in S is relatively easy, but computations for general 
number of cuts n require more systematic approaches. One such approach is to evaluate period 
integrals systematically; ref. [73] established a methodology, computing JF for general n up 
to S 3 terms. Another approach is to use the relation to matrix models. One can evaluate the 
matrix integrals directly [80] or by a more sophisticated diagrammatic technique [72]. This 
matrix model approach turns out to be rather efficient in actual computations and indeed, in 
section |D~T} we will compute JF up to S 5 terms. Yet another approach is to use the relation to 
the Whitham hierarchy [81]. For other work on computations of JF for general n, see [82-84]. 

We believe that the result of this appendix has various practical applications, including 
search for nonsupersymmetric vacua in M = 1 gauge theories. 

D.l Matrix Model 

By the Dijkgraaf-Vafa relation [63-65], the prepotential J-'q(S), S = (S 1 ,...,S n ), of the 
geometry fID.lj) is related to the free energy of the associated U(N) matrix model, 



Z = e - Fmm{9s ' N) = J d N \ exp 



— trW($) 

9s 



;d.4) 



where $ is an iV x N matrixjfj This matrix integral is performed around the vacuum where 
N % eigenvalues of <3> sit at aj. If we replace N % in F mm (g s , N) by S l by the relation 

9sN l = S\ (D.5) 

then the free energy organize itself into a genus ('t Hooft) expansion. Namely, 

F mm (g s , -) = f> s 2 *- 2 ^(S). (D.6) 

V 9sJ 9=0 



31 Here, the argument "TV" in F mm (g s ,N) denotes (N 1 , . . . ,N n ) collectively, not to be confused with the 
rank N — ^\ TV 1 of the matrix $. 
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As reviewed in section |5l one can evaluate the matrix integral (ID.4[) by perturbation theory 
using diagrams [72], as far as the perturbative part of F mm is concerned. However, this quickly 
gets out of hand, particularly because for general n one can have p-point interaction vertices 
with arbitrarily large p, which makes the number of diagrams explode. Namely, if we expand 
$ around the critical point x = 0,%, each coefficient g^ v in the expansion 

2 "+ 1 

W(oi + x) = W{on) + ^-x 2 + Y 9 -^x p (D.7) 

2 U p 



gives a p-vertex interaction, and p can be arbitrarily large for general n. Here 

1 
(P-1) 



m 2 = W"(a t ), g iip = —L^W<*\ai) (D.8) 



and W^ is the pth derivative. 

A more efficient method amenable to computer was proposed in [85-87], and here we 
generalize it to the case with an arbitrary number of cuts n. First note that the perturbative 
part of the matrix model free energy F mm can be written as an expansion in the coupling 
constant g s as: 

F mm>pert {g s ,N) = J2g k s f k (N). (D.9) 

fe=i 

Here, the order k amplitude fk{N) is a polynomial of degree k + 2 in iV J, s, which in turn has 
a genus expansion as follows: 

m 

f k (N) = J2 4":t 29+ / 41 • • • N lk -^\ (D.10) 

9=0 

where the coefficients A\? are totally symmetric in ii, «2, • - • , and [x] is the integer part 
of x. For a given finite k, the number of coefficients A\ '£ ^ n fk{N) is of course finite. 
Therefore, if we compute fk(N) for some small values of {N 1 } by computer, we can determine 
the coefficients A\fJ . Furthermore, there is symmetry under exchange of eigenvalues; for 
example, if we know A^fg, we can obtain A^u by ^ ne manipulation: 

(oil, mi, 9i, P ) <-> («2, m 2 , g 2 , p ), (a 3 , m 3 , g 3iP ) -^ (a 4 , m 4 , g AjP ). 

This symmetry significantly reduces the number of "data points" {N 1 }, for which we should 
evaluate the matrix integral in order to determine fk(N). In particular, this means that, if 
one knows fk(N) for n = k + 2 cuts, then one can determine fk{N) for arbitrary number of 
cuts n by symmetry. 
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For actually evaluating matrix integrals, it is convenient to go to the eigenvalue basis [88]: 



-Fmm(g,N) 



d N X 



N 



II(Aa - x b y 



.a<b 



exp 



1 N 

ys a=l 



(D.ll) 



where the Van der Monde determinant is from the change of variables [88] . We would like to 
compute this perturbatively around the vacuum where N* of the eigenvalues A a 's are equal 
to ctj, where i = 1, . . . , n. So, let us divide A a 's into n groups and expand around a, as: 



Xia = ca + Hia- i = l,...,n, a = l,...,N\ 
Then the matrix integral (ID.llj) is, up to a multiplicative constant, 



(D.12) 



d j y 



n 7V ! 



nn^ - ^)' 



i=l a<b 



n N i Ni 

I I II ^ia ~ Vjb + <*i 
i<j a=l 6=1 



x exp 



n N i i 2 

m i „2 i V^ 9i,p p 

(i ^ ^ \ 2 ia 2-^ v in 

ys i=\ o=l \ p=3 F 



n N* 

rEE 



(D.13) 



where we used the expansion ( 1D.7|) and 



aij = «j — aj. Given {N 1 }, we can evaluate this using 
computer by power expansion in g s which, following the procedure sketched above, allows us 
to determine fk{N) order by order. 



D.2 Result 

By setting Aj = mf, the first order result (0(N 3 )) is 



/i(*)=E( 



9iA 2 9 



U 



2A 2 3A? 



m 3 + J2 



2^,3 1 



A 2 o--- Ao- 2 . A-o- 2 . 



4£ 



i<i<fc v A i%'^ A ^%' A fc«H«ifc, 



(N l yN ] 



N l N J N k 



+ £ 



#i,4 0i,3 



N l 



(D.14) 



,4 A 2 6Af 

The coupling constants ^ iP can be expressed in terms of a 1 using (1D.8J) . The terms cubic in 
N l are planar amplitude, while the ones linear in N l are genus 1 (torus) amplitude. This 
agrees with the known result [73], if we use identities 

9%A 1 v^ 1 9i,3 ( 1 1 v^ M 

Jn+l-7- > , T"2 — <?n+l I 2 A "l A~ 2^ 

A- ^— » n/.-.-n/.-j. ni : : l\. \ /-vr./V- /-v-./V- ^— » ru-, 

(D.15) 



A 2 " ^ 






O^ijO^ik 



<XiAi 



afjAi aijAi ^. a ik 



A? ^ +1 , ^ a|A, ^ ^ oyo, J 
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upon using which flD.14[) becomes, after setting g n+1 = 1, 

l ^ Jjkl l J J ]<k J ' 

+ W— — + -A_ v —) nvw 

f-<\A t a^ A^% ottsbi^ctoP 



+ 4 y (—^ + — - + — — ) N l WN k 

i<~<k ^ AiCXijaki AjCHjkaij l\ k a ki a jk J 






j>k^i 



The second order result (0(iV 4 )) is much more lengthy. Let us write it as: 

f 2 (N) = J2 «««^ + Yl a m N i N J + E a ^ N i N I 
i i^ij i<j 

+ J2 a m N ? N i N i+ Yl a ijk iNiN 3 N k Ni 

i,j,k i<j<k<l 

j<k 

+ yb u N 2 + Y j b ij N i N r 



i<j 



Then the coefficients ana, etc. are: 



%,6 _ ^9i,59i,3 _ ^9i,A Q 9iA9i,3 _ 8gt,3 
«--+ 6A 3 A 4 g A 4+ A5 3 A6> 



'..3 



%,5 __ 12^,4^,3 _ 2g iA 8^3 4g 2 3 8ffj, 3 _ 8g jj3 _ Ag ifi 

""" " + a l3 A 3 ~ a l3 Aj a|A? + a,,Af + c^A 4 + 3oJA? Sc^A 3 a|A|A, 
1 4 4 



+ ^^T + 



a 4 A? a 4 A 2 a 4 AA-' 
6^,4 6ffj [4 8fl? 3 8g 2 3 2g ij3 g jt3 



a 2 A 3 a 2 A 3 a 2 A 4 a 2 A 4 a 2 A 2 A 2 

8&,3 8ffj, 3 2ffj |3 2^3 5_ 5 11 

a 3 A 3 a 3 A 3 . a? A 2 A ■ a 3 A-A 2 a 4 A 2 o; 4 A 2 «*.A-A- ; 
12^. 16c/| 



O'iijk a3 a 4 



8ft, 3 8ft i3 8<7j t 3 8ft ;3 4ft ;3 4ft ;3 

a%a ik A 3 a 2 k a l3 A 3 a 2 j( x jk A 3 a 2 k a jk A 3 a^a^AfAj a 2 k a jk A 2 A k 



+ —, T7T + 



oiija ik a 2 k AjA k a^a^A 2 a^a^AiAj a 3 k a jk A 2 a^ayAfA 
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aJaifcA? afjOtjkAiAj af k a {j A 2 af k a jk AiA k 
4 2 4 

+ a 2 -a 2 ,A 2 ~ a 2 a 2 A 2 + a 2 , a 2 , A 2 ' 

16^,3 16^.3 16# fe , 3 16#, : 

Qijfc/ — + 



oiija ik oiiiA\ aijOLjkOijiA^ a ik a jk a kl A 3 k auajia k iAf 



+ 



ctfjOiiiaijkAiAj afjUikajiAiAj ' a 2 k aua jk AiA k a 2 k aija k iAiA h 

8 
+ ^ ^-^ + ^ ^-^ + ^ ;— ^ + 



a^anzajiAiAi a^a^akiAiAi a 2 k a ik ajiAjA k a 2 k aija k iAjA k 



+ 



a^aijakiAjAi ajpua^AjAi a 2 kl aua jk A k A t a 2 kl a ik a ji A k A t 
8 8 8 8 

a^onkOiiiA 2 a^ajkajiA 2 a^a^auA 2 a 2 k a jk a kl A 2 
8 8 8 8 

+ ^ 7^7 + 



a^a^aikA 2 auajiakiAf a^a^ajiA 2 a 2 k a ik a k iA 2 k 



+ —, ^ + 



a^aijCXjkA 2 a^a^a^Af a kl a ik a jk A 2 k a 2 kl a a ajiA 

"~ + 3A 3 A^ + 2Af 8A? 3Af 

, 2 &,5 2ff j|5 ^9i,A9i,3 4:9j,39jA AM 9jA 

"ii~r a-3 ao » a T 



ayA? OijAJ a^Af otij-Aj a|Af a^A 3 
2 ft ? ,3 2 ^,3 9i, 3 9).3 2g i>3 2^-3 

ac; a k ' oa/1 oa/1 



c^Af ayA} aJA* a?. A} Sc^A 3 So^A 3 
2 11 

"4A^ + 24Af + 24Af (D ' 18) 

The third order result (0(iV 5 )) is too lengthy to be included here. The interested reader 
can find the result in the Mathematica file included in the source file for the current paper at 
arXiv.org. 

It is easy to identify the matrix model diagrams corresponding to each term in the above 



result. For example, the first term in ana in the second order result (ID. 181) . -^ comes from 
the following two planar diagrams: 
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Note that the definitions (1D.4|) means that the coupling constants g^ p enter the free energy 



Fmm with a sign as follows: 

-IB-**)- ( D - 19 ) 



i,p 



This is in addition to the signs coming from the fermionic ghosts. 
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